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Abstract 



Starting from general considerations concerning phase transitions and the many- 
body problem, a pedagogical introduction to the Mermin- Wagner theorem is 
given. Particular attention is paid to how the Mermin- Wagner theorem fits in 
with more general methods in many-particle physics. An attempt at an assess- 
ment of its validity for approximative methods is made by contrasting the results 
obtained within Onsager reaction-field theory and those obtained by the well- 
known Tyablikov method. 
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1. 



Introduction 



1.1 



Outline 



In this paper, we shall investigate the role of spatial dimensionality in the context 
of many-body theories for magnetic long-range order. After a rather general in- 
troduction to the problem (this chapter), and an exposition of the basic notions 
of phase transitions and spontaneous symmetry-breaking (section 2.1), the Bo- 
goliubov inequahty is discussed (section 2.2). This is a rigorous relation between 
observables of a physical system and the system's Hamiltonian, and it shall be 
used at length in Chapter III in the course of an extended proof of the absence 
of a magnetic phase transition in films described by the Heisenberg, Hubbard, s-f 
(Kondo-lattice) and Periodic Anderson model. Chapter IV investigates the con- 
nection of the Bogoliubov inequality and related correlation inequalities with the 
Green's function technique widely used in many-body physics. Finally (Chapter 
V), two examples of approximative theories are discussed to illustrate the theo- 
retical findings of the previous chapters. 

1.2. The many-body problem 

Loosely speaking, the success of physics in describing properties and processes of 
the physical world is, to a large part, the result of taking a reductionist viewpoint, 
in the sense that macroscopic phenomena are regarded as being fully determined, 
at least in principle, by underlying microscopic processes which behave according 
to some (more or less well-understood) fundamental law of nature. In the case 
of electrons in solids, one would expect the non-relativistic Schrodinger equation 
(plus the Pauli principle) to provide an accurate description of the behaviour of the 
physical system, at least in the low-energy regime. For non-interacting electrons 
in a completely periodic potential of ions located at the lattice sites of a perfectly 
ordered crystal, this can indeed be accomplished. However, it turns out that, 
for the many-body problem of strongly interacting electrons, such an approach is 
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highly impracticable, and attempting to calculate macroscopic properties from a 
set of 10^'^ one-particle equations simultaneously, would be fighting a losing battle. 

One solution to this problem is to propose reduced many-body models which 
retain the basic features of the microscopic dynamics, e.g. Coulomb repulsion 
or the Pauli exclusion principle, but at the same time leave out less important 
features, e.g. by neglecting interactions between electrons and phonons or by 
limiting Coulomb repulsion between electrons to those at the same lattice site. 
In most instances, such models still cannot be solved exactly; methodologically, a 
"computational gap"Q remains which one can only hope to close, to some extent, 
by use of approximative methods, whose validity must be constantly reassessed. 

1.2.1. Summary of many-body Hamiltonians 

The main focus of this paper is on magnetic phenomena. It can be rigorously 
shown (Bohr-van Leeuwen theorem) that the existence of magnetic phenomena 
(diamagnetism, paramagnetism, and collective magnetism) of an ensemble of par- 
ticles in a crystal is inconsistent with classical statistical mechanics.^ This striking 
fact indicates that a satisfactory description of spontaneous magnetic order will 
have to include quantum effects from the very beginning. In the following, we 
shall discuss almost exclusively quantum many-body models. Nevertheless, clas- 
sical models such as the Ising modelQ defined by the Hamilton function 



where Si is a classical (c-number) "spin" at lattice site i, have played, and continue 
to play, an important role in the conceptual development of the theory of phase 
transitions. They serve as a test bed for new mathematical approaches, such as 
the renormalization group approach, and as limiting cases of more sophisticated 
models. 

The Ising model, for example, can be derived from the more general Heisen- 
berg model as defined by the Hamiltonian 



^For a concise discussion of this fact see [Q. 
■^Also, and more appropriately, known as Lenz-Ising model, to give credit to the true initiator 
of the model, W. Lenz (for a broad account of the history of the subject, see H]). 







^This expression was coined by j^]. 



by singling out one of the components (x, or z). The Heisenberg model is a 
quantum many-body model and the spin Si must now be interpreted as a quantum- 
mechanical operator. From the x and y-components of Si one constructs spin 
raising and lowering operators 

Sf = S!±tSf. (1.3) 

The operators Si and 5*,^ obey the commutation relations for angular momentum 
operators so that, among others, the following identities hold: 

[S^,Sf]^ = ±h5,,St (1.4) 
[St^Sj]^ = 2h5,,St (1.5) 



and 
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Si ' Sj — — (^S^ Sj -\- Sj^ Sj ) -|- Si Sj . (1.6) 

Note that the indices i and j are taken to identify uniquely the lattice sites of 
spins Si and Sj. Instead of numbering all sites consecutively, one could identify a 
lattice site by more than one index, e.g. in film systems where it will turn out to 
be useful to identify lattice sites by numbering layers a,p... and positions i,j, ... 
within the layers separately: i i— >■ {i,a). The Kronecker symbol Sij in ( |1.4| ) and 
( |1.5|) must then be replaced by the product dijd^p. 

At times it will be helpful to make use of Fourier transformation into the 
reciprocal lattice: 

i 

with the inverse transformation 

g{+,-,-,y,z) = 1 J] ^-^msi+^-^-^y^^){k). (1.8) 

k 

The Heisenberg model was proposed some 75 years ago [|], in order to de- 
scribe the behaviour of magnetic insulators such as EuO (ferromagnetic), EuTe 
(antiferromagnetic) and EuSe (ferrimagnetic). Whether a Heisenberg- type mag- 
net will order ferromagnetically or antiferromagnetically depends on the sign of 
the coupling constants Jjj, also known as exchange integrals. When Jij > 



{Jij < 0) for neighbouring spins Si, Sj, ferromagnetic (antiferromagnetic) order 
will be favoured, as this will minimize their contribution to the total energy of the 
system. Different exchange mechanisms exist; the feature they have in common is 
that by virtue of quantum mechanical symmetry requirements, electrostatic cor- 
relations lead to an effective spin-spin interaction. The numerical values of Jij 
depend on details of the substance and the exchange mechanism involved. 

The Hubbard model, in contrast to the Heisenberg model, describes the 
case of itinerant electrons in narrow energy bands, as found in the prominent 
ferromagnetic metals Co, Ni, Fe. It was proposed in 1963 by Hubbard ^ (and 
simultaneously by Gutzwiller and Kanamori [§) and takes into account that 
in these systems the same group of electrons is responsible for magnetic effects 
and electrical conduction. Further empirical findings that are taken into account 
in the Hubbard model are, firstly, that narrow 3(i-bands are responsible for the 
appearance of collective magnetism; secondly, the screening of the Coulomb inter- 
action between d-electrons through electrons of broad s- and p-bands; and thirdly, 
the importance of the lattice structure (for a detailed derivation see The 
individual points can be discussed by looking at the structure of the Hubbard 
Hamiltonian in its simplest (one-band) form 

ija i(T 

where and cja are creation and annihilation operators for an electron with 
spin cr at lattice sites i and j, respectively; and ni±„ = cf^^Ci±a- is the number 
operator which tests whether or not an electron with spin icx is present at lattice 
site Z.0 The Hubbard Hamiltonian consists of a kinetic part and an interaction 
part reflecting the Coulomb repulsion between electrons which becomes relevant 
only if two electrons happen to be at the same lattice site (which, because of the 
Pauli principle, must then have different spin quantum numbers). The hopping 
integrals Tij, i.e. the probability amplitude for an electron with spin a at site j 
to move to another site i, are related to the one-particle Bloch energy dispersion 
e{k) by Fourier transformation: 

Ti, = l^£(fc)e^^-(^»-^^) (1.10) 

k 

«i 

^For a detailed account of second quantization in many-body physics, see [pj. 



Spin-like operators, i.e. pseudo-spins a, are built from the electron creation 
and annihilation operators, in the following way: 

at = hc^c,^ (1.12) 

= ^c+^^T (1-13) 

Fourier transforms into reciprocal lattice are defined in exactly the same way 
as in ( lOD and ([TSl ), e.g. 

a+{k) = hJ2e''-'^'ctc,^- (1-15) 

i 

In the s-f (Kondo-lattice) model and the Periodic Anderson Model (PAM), 
the situation is somewhat different from the Heisenberg and Hubbard models, as 
one must now distinguish between two spin systems. Both models deal with a 
situation that can be found in the rare earths, where the electrons of the inner 4/ 
(5/) shells, which are incompletely filled and which correspond to strictly localized 
states, are screened by the broader 6s (7s) states. Thus, the magnetic moments 
associated with the 4/ (5/) electrons couple according to Hund's rules and form 
a local magnetic moment. The s-f model describes the coupling of such a local 
/-spin S to the electrons of the conduction band: 

ij(T i(T 

where z^^r = ±1. In the PAM the Hamiltonian includes the hopping Hg of the 

electrons of the conduction band (as in the s-f term Hj defining the 

energy level of the /-electrons, the Coulomb repulsion Hff between electrons in 
the strongly locahzed /-states, and the hybridization Hhyh of the /-level with the 
conduction band: 



H = H, + Hf + Hff + Hhyb (1.16) 

U 



2. Phase Transitions and the 
BoGOLiuBOv Inequality 



Among the general features of physical systems at a macroscopic scale, the display 
of different phases plays a fundamental role. By this we mean that the overall 
characteristics of a physical system, for example a large bulk volume of a certain 
substance, can be described by a limited number of system-specific parameters 
(in many cases of interest by a single one) which themselves depend on a set 
of certain external conditions, usually an external field, and the temperature. 
When a system changes from one phase to another, e.g. as a result of increasing 
temperature, this transition will often be marked by a rapid change in one of 
the system's characteristic quantities, such as the density in the liquid-vapour 
transition. From early on, it has been considered a challenging task to describe the 
experimental findings through theoretical models. On a phenomenological level as 
well as by starting from first principles, substantial progress in the description of 
classical systems has been achieved since the pioneering work by, for example, van 
der Waals and Weifi in the liquid- vapour and the ferromagnetic phase transition. 
With the advent of quantum theory and progress in the understanding of the laws 
governing the microscopic behaviour of physical systems, it became possible to 
"design" mathematical models in such a way as to reflect certain aspects of the 
microscopic dynamics of real systems while omitting others. This has led to a 
remarkable wealth of models, some of which have been described in some detail 
in the previous chapter. 

2.1. Phase transitions and spontaneous symmetry-breaking 

In the theory of condensed many-body systems one is often interested in de- 
scribing phases that fail to exhibit certain symmetry properties of the underlying 
Hamiltonians. This is, in fact, a common feature and one can classify such systems 
according to their symmetry properties. Thus, crystals, by their very lattice struc- 
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ture, break the translational symmetry encountered in the continuum description 
of fluids; ferromagnets, in addition to the spatial symmetry-breaking due to their 
crystal structure, are not invariant towards rotations in spin-space, even though 
the underlying Hamiltonians describing the system may well be. Less obvious 
types of symmetry-breaking occur in other quantum systems, such as superfluids 
and superconductors, where gauge invariance is broken. 



2.1.1. Bogoliubov quasi-averages 



Bogoliubov |T3)[|ni has devised a method to describe the occurrence of sponta- 
neous symmetry-breaking in terms of quasi-averages^ Normally, for systems in 
statistical equilibrium, we have defined the average value of A as the trace over the 
equilibrium density (or statistical) operator p = exp(— /37i) times the observable 
A, so that for the infinite system {V —>■ oo) one has 



{A) = lim tr (pA) 

V^oo 

= lim tr {e'^'^A) , (2.1) 



where Ti. is the grand-canonical Hamiltonian, Ti = H — pN {N: number opera- 
tor). However, it turns out that under certain conditions such averages may be 
unstable with respect to an infinitesimal perturbation of the Hamiltonian. If a 
corresponding additive contribution H^, = vH' of the order of v is added (where 
1/ is a small positive number which will eventually be taken to zero: 0), i.e. 

n, = H + H,- pN, (2.2) 

one can define the quasi- average of A in the following way: 

(A)^ = lim lim tr(e"^^'^A). (2.3) 



The average ( p.l|) and the quasi-average ( |2.3| ) need not coincide, since the two 



limits in ( |2.3| ) may fail to commute within some parameter region (i.e. for some 
combination of p and /?). 

Quasi-averages are appropriate for cases in which spontaneous symmetry- 
breaking occurs, as can be shown by a simple argument. Suppose the Hamiltonian 



^This section follows [10 and [h2 



Ti displays a continuous symmetry S, i.e. it commutes with the generators of 
the corresponding symmetry group, 

[H,ry_=0. (2.4) 

If some operator B is not invariant under the transformations of iS, 

[i?,r^]_ = cvo, (2.5) 

the (normal) average of the commutator vanishes, 

{C') = 0, (2.6) 



as can be readily seen from eqn. (|2.1| ) by use of cyclic invariance of the trace. 
In those instances, however, where the perturbative part Hi, of ( p. 21) does not 
commute with F^, this will give a non-vanishing quasi-average: 

= limtr(e-^^^[5,ry_) 

^ 0. (2.7) 

Thus, even though one might naively expect the quasi-average to coincide with 
the "normal" average in the limit — > 0, quite generally this will not be the case. 
Note that the quasi-average depends on the nature of the perturbation added to 
the "original" Hamiltonian. 

In the following, we will make extensive use of the concept of quasi-averages 
as in ( p.3|) .p| As a first example, let us examine the Heisenberg model 

H = -Y,J^J{S,■S,>|, (2.8) 

ij 

which is invariant under the continuous rotation group generated by the total 
spin vector S = J2i because of H, S = 0. Thus, we can take S as the 

operator B, and from ( p.6|) one finds ^ [5'", S*^] _^ = 0, where S"", S"^ are the 

components of the total spin vector S. Together with the commutation relations 
for spin operators, e.g. [S^,Sy]_ = ihS^, one sees that the conventional average 

^The symbol (...) will be used for averages and quasi-averages alike, the distinction between 
the two being obvious from context in all cases under consideration. 



of the magnetization vanishes. This is just a manifestation of the fact, that on the 
macroscopic level, for an ideal, infinitely extended system, there is no preferred 
direction in space. One can say this is a situation of "degeneracy" with respect 
to spatial orientation. Adding an external field, Bq = BqCz along the z-axis 
for example, lifts this degeneracy and via the contribution H}, ~ BoM{T, Bq) to 
the Hamiltonian (M: magnetization) one can, thus, construct appropriate quasi- 
averages for arbitrary operators according to equation (|2.3| ). 

2.1.2. Order parameters 

In the theory of phase transitions, the first step is to identify a quantity whose 
(quasi-) average vanishes on one side of the transition, but takes a finite value on 
the other side. This quantity is called the order parameter. In a continuous phase 
transition, the order parameter may gradually evolve from zero at the critical point 
to a finite value on one side (usually the low-temperature side) of the transition. 
For different kinds of phases, different order parameters must be chosen. From 
a phenomenological point of view, one must consider each physical system anew. 
We have already mentioned the liquid- vapour and the ferromagnetic-paramagnetic 
transition as prototypes for phase transitions. In the former, the obvious choice 
for the order parameter would be the difference between the mean densities, i.e. 
P — Pvapour', in the latter the relevant order parameter is the magnetization M. 
Within a given many-body model, the magnetization can be defined in microscopic 
terms, as we shall see shortly. 

Sometimes, as in the transition to the superconducting state, it may even 
be possible to characterize the same type of phase transition by use of different 
order parameters. According to the standard theories of superconductivity, at low 
temperatures electrons with opposite spins form Cooper pairs. Thus, a possible 
order parameter would be the average probability amplitude to find a Cooper pair 
at a given lattice site in the crystal. Alternatively, one could characterize the phase 
transition through the gap parameter whose modulus is the difference in energy 
per electron of the Cooper pair condensate and the energy at the Fermi level. 
We shall briefly return to the problem of superconductivity when we consider the 
possibility of pairing at finite temperatures in film systems. 



2.2. Bogoliubov inequality 



2.2.1. Proof of the Bogoliubov inequality 



The Bogoliubov inequahty is a rigorous relation between two essentially arbitrary 
operators A and B and a valid Hamiltonian if of a physical system. In its original 
form, proposed in [|ll|], it is given by 

\{[C, ^]-> r < f ( ^1 +) ( [C^ [H, C]^]_) , (2.9) 

where (3 = l/ksT is the inverse temperature and (....) denotes the thermody- 
namic expectation value. A and B do not necessarily have an obvious physical 
interpretation from the very beginning, so the physical significance of ( p.9|) will 
depend on the suitable choice for the operators involved. 

The inequality can be proved by introducing a scalar product which is based 
on the energy eigenvalues and the corresponding energy eigenstates \n) of the 
Hamiltonian H : 



■'nm 



{n\m) = 6„ 
H\n) = En\n) (2.10) 
^ {n \H\ n) = En 

with the completeness relation (1: identity operator) 

l = J2\n){n\. (2.11) 

n 

By further introducing the statistical weight 

eM-PEn) 

" tr{expi-PH)y ^ ' 

the Bogoliubov scalar product 53 of two arbitrary operators A, B is defined as 

^{A-B)= (2-13) 

n,m 

En^Era 

Note that the sum is restricted to terms with En 7^ Em only. Since 03 is positive 
semidefinite (see or appendix A), the Schwarz inequality holds: 

|Q3(A; B)\^ < Q3(A; A)Q3(5; B). (2.14) 



With the specific choice 

B = [C+,H]_ 
the mixed product 03 (A; B) on the LHS of ( glil) is 

^{A;B) = ^ {n\A+\m)(m [C+,H] 

n,m 



(2.15) 



n 



E - E 



n,m 



(2.16) 

by use of the completeness relation ( ^.11[ ) and the definition of 

The norm ^{B-B) = ||S|r > can be calculated from eqn. ( p.l6|) by spe- 
cializing A = B = H]_ , thus arriving at 

^{B;B) = ([C^,[H,C]_]_)>0. 



(2.17) 

For the remaining norm Q5(/l; A) = \ \A\f we start from ( [^.13D and note that 
the quotient {Wm — Wn)/{En — Em) > {Em 7^ En) can be bounded from above: 



En Em 



-I3E„ 



+ e 



) 



-PEm p — fiEn 



tr{exp{-(3H)) En - Em e"/^^™ + e'/^^- 



{Wm + Wn) 
En Em 



tanh [-{En- Em) < - {Wm + Wn) , (2.1 



where the condition Em 7^ En has again been used. It is now straightforward to 
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show that the norm \ \A\\ is bounded from above: 

^{A-A) < ^ {n\A+\m){m\A\n){wm + Wm) 

n,m 

Eny^Em 

- § l^'^ I ("^ 1^1 ^) {Wm + Wm) 



IE 



Wm. m 



[AA-l 



m 



(2.19) 



This is all that is needed for deriving the Bogoliubov inequality; inserting ( |2.16| ), 
( p.lTp , and (|2.19|) into the Schwarz inequality (|2.14|) immediately recovers (after 



minor rearrangement of the LHS) the Bogoliubov inequality in the form proposed 
in eqn. (|2.9| ). From the derivation it is noteworthy that the two factors on the 
RHS each are upper bounds to a norm and, thus, greater than or equal to zero. 
In particular, if, for example, the double commutator depends on some parameter 
k, one will always have 

[[C, HUC^]_) {k) + ( [[C, HUC^]_) {k') > ( [[C, i/]-, C+] _) {k). (2.20) 

In later calculations a slightly modified version of the Bogoliubov inequality 
will turn out to be useful. Dividing both sides of the Bogoliubov inequality ( |2.9| ) 
by the double commutator and summing over all wave vectors k associated with 
the first Brillouin zone in the reciprocal lattice, one arrives at 

y — < ^ y / u A+l \(^). (2.21) 

2.2.2. Applications of the Bogoliubov inequality 

Hohenberg [|l^ was the first to note that the Bogoliubov inequality could be used 
to exclude phase transitions, showing that there could be no finite-temperature 
phase transition in one- and two-dimensional superfluid systems. At roughly the 
same time Mermin and Wagner, following Hohenberg, considered the case of spon- 
taneous magnetization in the Heisenberg model. Incidentally, the Mermin-Wagner 
theorem, as it became known later, was published shortly before Hohenberg's 
proof, and it has since become a landmark paper in the theory of phase transitions 
[|T^ .0 A variety of extensions of the theorem as to the possibility of spontaneous 
magnetization have been given (see introductory remark in the next section) for a 
number of many-body models. Here, we shall briefly mention further applications 
to other order parameters.^ 

Mermin applied the Bogoliubov inequality to the possibility of crystalline 



order in two dimensions, confirming earlier suggestions by Peierls and Landau 



^Note that Mermin and Wagner never claimed priority, instead explicitly mentioning Hohen- 
berg's suggestion in their paper. 

survey of all papers concerning extensions of the Mermin-Wagner theorem is forthcoming 
in a review article by this author, see ref. pi]. 



that there could be no two-dimensional crystalline ordered state. The proof is 
carried out in detail for classical crystals, giving only an outline of the quantum 



case. For the latter case, the derivation was given by Fernandez [16 



Kishore and Sherrington considered a quite general Hamiltonian of elec- 
trons interacting non-relativistically among themselves, and with spatially ordered 
or disordered scatterers, excluding spontaneous low-dimensional magnetic order. 
The restriction to non-relativistic interactions means that spin-orbit effects are 
excluded from the problem so that the results keep the backdoor open for possible 
Ising-type, or other suitable anisotropics. 

The problem of partially restricted geometries was considered by Chester et al. 



18[| for three-dimensional Bose systems of finite cross section or thickness. The 
results of the generically one or two-dimensional case are, in principle, reproduced. 
Their main line of argument rests on the operators in the Bogoliubov inequality 
being defined as Fourier transforms in a restricted space DxDxL or DxLxL, 
respectively, where D remains finite in the thermodynamic limit while L goes to 
infinity. This way, the operators lack the physical meaning of the usual real-space 
or fc-space operators, which means that an intuitive physical interpretation is not 
readily obvious. 



3. Absence of Finite-Temperature 
Long- Range Order in Partially 
Restricted Geometries 



Following the well-known 1966 paper by Mermin and Wagner [M|, in which the 
authors prove the absence of spontaneous magnetization in the one- and two- 
dimensional isotropic Heisenberg ferromagnet (and the two-sublattice antiferro- 
magnet), much research has been carried out in order to extend the statement 
to other spin systems. Wegner considered a model describing a system with lo- 
cally interacting itinerant electrons |19|. Walker and Ruijgrok discussed a band 
model for interacting electrons in a metal Ghosh more specifically, re- 
covered the Mermin- Wagner theorem for the Hubbard model. A proof for the 
s-d interaction model is given by van den Bergh and Vertogen A paper by 
Robaszkiewicz and Micnas p3[ extends the Mermin- Wagner result to a general 
model with localized and itinerant electrons, covering the modified Zener model, 
the extended Hubbard model and s-d models as particular cases. In addition to 
extending the Mermin- Wagner theorem to several microscopic many-body models, 
generalizations to more complicated geometries were also found. Baryakhtar and 
Yablonskii ||2^, for example, proved the Mermin- Wagner theorem for systems with 
an arbitrary number of magnetic sublattices, also excluding non-coUinear magnetic 
order when an external field is applied. Thorpe p5| considers the case of ferro- 
magnetism in phenomenological models with double and higher-order exchange 
terms. Similar results were obtained in the multi-sublattice case by Krzemihski 
26[| . Recently, Matayoshi and Matayoshi have discussed models with n-th 



nearest neighbour exchange interactions and tried to extend the Mermin- Wagner 
theorem to anisotropic exchange interactions; their proof, however, requires ex- 
tremely special conditions on the parameters of the model which appear to be of 
no physical significance. Also, attempts to extend the Mermin- Wagner theorem to 
special cases in three dimensions remain somewhat inconclusive p8|. More recent 
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papers will be discussed later in this chapter.^] 

The general procedure is essentially the same for most of the papers above: The 
Bogoliubov inequality is used with suitable operators defined in such a way as to 
give an upper bound for the desired order parameter, i.e. the (bulk or sublattice) 
magnetization. In general, one arrives at an upper bound that, in the limit Bq — > 
0, behaves in one or two dimensions as ~ {MBqY'^ or ~ l/(lnM5o)^/^ thus 
excluding a finite value of the order parameter (M denotes the bulk magnetization 
and Bq is the external magnetic field). In three dimensions, no such behaviour 
is found, i.e. the usual line of reasoning fails and spontaneous magnetization 
cannot be ruled out. Within the usual scheme following Mermin and Wagner, the 
dimensionality of a system enters into the calculation only through the volume 
element when integrating the final version of the Bogoliubov inequality over k- 
space. Thus, a microscopic picture of the transition from two to three dimensions 
does not emerge. 



3.1. Films with finite thickness 

As a result of the rapid progress of thin-film technology in recent years, one can 
now prepare and study systems with restricted geometries, such as films, in great 
detail. Thus, important parameters, such as the Curie temperature, magnetization 
and susceptibility, can be measured and discussed as functions of the number of 
layers d in a magnetic film. Experiments indicate |^| that, for real systems, the 



transition from 2D to 3D behaviour of, for example, the critical exponent (3, occurs 
within a narrow crossover region of d. The critical temperature also shows a 
strong d-dependence, with Tc{d) quickly approaching the bulk value Tc(oo) as d 
increases 
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With these observations taken into account, it appears desirable to improve 
one's theoretical understanding of how the transition from the monolayer to the 
bulk occurs. By referring only to the truly microscopic properties of the respective 
many-body model, it should be possible to study in detail the transition from 
two to three dimensions. Work in this direction has been done by studying film 
systems with symmetry-breaking contributions to the Hamiltonian, thus making 
possible the study of the Curie temperature as a function of d and the anisotropy 



parameter (for analytic approaches to Heisenberg films see for an extensive 
review of the magnetic properties of thin itinerant-electron films see [^). 

iRefs. [01, l3l, i7|, m. 



At a more fundamental level, the validity of exact results, such as the Mermin- 
Wagner theorem in two dimensions, to thin films may be tested. In this chapter, a 
proof of the Mermin- Wagner theorem for film systems within the main many-body 
models (Heisenberg, Hubbard, s-f (Kondo-lattice) model, and Periodic Anderson 
Model) will be given. 

3.1.1. Film Hamiltonians 

For reasons of simplicity, we shall restrict our attention to film systems composed 
of d identical layers stacked on top of each other. The calculations can easily be 
generalized to account for more complicated geometries. One can think of the film 
geometry as consisting of a two-dimensional Bravais lattice, i.e. the first layer, 
with N lattice sites and a ci-atomic basis that corresponds to the d layers being 
stacked up. Each lattice vector then decomposes into 

Ria = Ri + ra, (3.1) 

— * 

where Ri is a vector of the Bravais lattice and fa is the basis vector pointing 
to the a-th layer. In all of the following calculations, Greek indices label layers 
and Roman indices refer to sites of the Bravais lattice. It should be noted that 
translational invariance can only be assumed within each layer. E.g., the notion 
of a reciprocal lattice only makes sense when referring to the two-dimensional 
Bravais lattice. A similar statement applies to Fourier transforms being thought 
of as connecting real-space quantities with those defined in wave-vector space. 

In the following we shall discuss the Mermin- Wagner theorem for the Heisen- 
berg model, the Hubbard model, s-f model (Kondo-lattice model), and the Peri- 
odic Anderson Model (PAM).0 For film systems one has to distinguish between 
Bravais lattice indices and layer indices for all site-dependent quantities, such as 
spin operators S^^' '^'^'^\ annihilation and creation operators Cj-^"* or coupling con- 
stants Jj"^ which depend on two lattice sites. In this notation, the Hamiltonian 
for Heisenberg films is given by 

H = — J^j^^i^S^Sji^ + SlaSjp) — e ^'Sia, (3-2) 

ija(3 ict 

where the term b ^ e~^^'^^S^^ is due to the interaction with an external magnetic 
field b = ^J^^-^p . The interaction with an external magnetic field leads to the 

^Parts of this proof have been pubhshed in jBSf . 



magnetization 



*^-]^^E^-*«««>4EM„. (3.3) 

ia a 

where the phase factor e~^^'^^ accounts for both ferromagnetic and antiferro- 
magnetic ordering, depending on the choice of K. We also assume the couphng 
constants J^^^ as satisfying some general properties, i.e. 

= J^- (3.4) 
Ju = 

and 

2 



1 ^ ^ , \Rm- Rpj 
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The last condition ( |3.5[ ) is very weak considering that the exchange integrals J 
will usually decay exponentially with distance. 

As an example for itinerant-electron systems, we shall discuss the Hubbard 
model 



ija/3cr iaa ia 

where T^"^ describes the hopping of an electron from lattice site j of the /3-th layer 
to site i of the a-th layer and U is the energy associated with having two electrons 
at the same lattice site. Again, a term corresponding to an external magnetic 
field is included. The z-component of the pseudo-spin, a^^, is defined according to 

-iaf3 



( |1.14| ). Similarly to the Heisenberg case, we require the hopping constants Tj" to 

rpf3o 

'ij ~ ji 



satisfy the isotropy conditions T^^f = T^^ as well as to converge upon summation 



over all lattice sites: 

(jin — Rh 



7/3 nk 



= q <oo. (3.7) 



In the s-f model one deals with two spin sub-systems, ISiaj and {aia}, the 

former consisting of locahzed /-electrons, the latter being associated with itinerant 
s-electrons: 



(3.8) 

A similar situation occurs in the Periodic Anderson Model, where one must 
distinguish between the electrons of the conduction band (described by fermionic 
creation and annihilation operators c-^^) and the /-electrons (whose creation and 
annihilation operators we shall denote by fl^l)'- 

ijaficr iaa ia 

+V (c-^/- + f^^c..c^a.) - b e-^^'-^^^ + ) (3-9) 

iatj ia 

= Hs + Hf + Hff + Hhyb + H,. (3.10) 

Here, a^_^ and ctj.^ are the 2;-components associated with the pseudo-spins of the 
conduction and the /-electrons, respectively: 

= 2 '-'^^^"T ~ ''^ial') ^ '^Uta^fia} " f^lfial)- (3-12) 

3.2. Absence of a magnetic phase transition in film systems 

3.2.1. Choice of operators and evaluation of the Bogoliubov inequality 

By inspecting eqn. ( p.9|) one can easily conclude that the choice of suitable oper- 
ators A and C is crucial; it determines whether the inequality will be physically 
meaningful or not. In film systems, long-range magnetic order within a given 
layer is conceivable where the bulk (or even sub-lattice) magnetization of the 



whole system vanishes, e.g. 



• •• nnnnnnnn... 

• •• nnnnnnnn... 

Excluding long-range magnetic order for every layer within a film would, there- 
fore, be a considerably stronger statement. The general idea of our proof is to 
use the Bogoliubov inequality to find an upper bound for the layer magnetization 
Mp, i.e. 

M^</(i?o,M), 

where / is a function that approaches zero as -Bq ^ and does not depend on 
any layer-specific quantities. 

This is best achieved by choosing C as a sum Ylia^---) spin operators and A 
as a spin operator associated with a specific layer (3. Hence, the numerator of the 
LHS of the Bogoliubov inequality will be a layer-dependent quantity, while the 
double commutator in the denominator will be summed over all lattice sites and 
thus be layer- independent. One may then expect to be able to replace the RHS 
of the inequality by a (layer- independent) upper bound. More specifically, for the 
Heisenberg model we set 

A^^)=S~{-k-K) (3.13) 

and 

I — * 

where the Fourier transforms Sp{k) of the layer-dependent operators are defined 
within the two-dimensional reciprocal lattice according to eqn. ( |1 . 7|) . Note that 



in fc-space we have (^S^{k)^ = S*^ as is obvious from the definition of the 

Fourier transform. The operators S in the context of the Heisenberg model are 
spin operators for which the usual commutation relations (|1.4| ) and (|1.5|) hold 
true, which, after Fourier transformation into fc-space, run 

S+ih),S0{h)] =2h5^pSl{k, + h) (3.15) 



and 



= ±hdapS^{h + k2). (3.16) 

In the Hubbard model a similar definition for the operators A and C is used. 
However, as mentioned above, the operators S^^^^'' are now built up from fermionic 
creation and annihilation operators, e.g. 

a-(fci) = hci^^^c,;^^^ (3.17) 

and 

4(4) = hcl^.c^,^^^. (3.18) 
The operators A and C are, then, 

Aa) = CT-i-k -K)^ f^C^_U_ji^^^C_^k-KM (3-19) 

and 

c ^Y.^t{k) = hY^ci^^c^m f= E^/^) ' (3-20) 

respectively. The commutation relations for spin operators may, in a purely formal 
sense, be used in this case as well. 

In the s-f model and the Periodic Anderson Model, one must not forget that one 
is dealing with two distinct spin sub-systems. In the former, one of the subsystems 
can be described by ordinary spin operators, while the other is associated with 
itinerant electrons; in the PAM both subsystems must be described in terms of 
fermionic creation and annihilation operators. In the s-f model, the two systems 
are independent from one another in the sense that spin operators and creation 
or annihilation operators commute: 

= (3.21) 
and, thus, 

= 0, (3.22) 
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{+,-,x,y,z) (+) 



(+) 



while in the PAM similar anticommutation relations hold between creation and 
annihilation operators of the conduction electrons and the /-electrons: 



c, 



(+) A+) 

iaa ' J ji3a' 



fciatr' fc2/3o-' 




0. 



With (|3.18|) ,( ^yT^ ) we define as operators A and C in the s-f model 



(3.23) 
(3.24) 



and 



/3 fS 

and in the Periodic Anderson Model 



A(^) = S-{-k -K) + a-{-k - K) 



^7) = <y-4-k-K) + af^{-k-K) 



(3.25) 



(3.26) 



(3.27) 
(3.28) 



3.2.2. Hamiltonian-independent quantities 



It is obvious from the structure of the Bogoliubov inequality ( p.21j ) that the nu- 
merators on both sides of the inequality are determined entirely by the choice of 
operators A and C as mentioned earlier. The Hamiltonian H of the many-body 
model enters the calculation only via the double commutator ([[C, C"'"]_) . 
For simplicity, we shall start with the Hamiltonian-independent expectation val- 
ues appearing in the Bogoliubov inequality, more specificially with the [C, 
commutator 



{ [S^ik) + a^ik), S-i-k - i?) + a-{~k - K) 

mnf3 



(3.29) 



for the s-f model. For reasons described above, S and a operators commute, so 
the commutator can be evaluated directly by use of the fundamental commutation 
relations. This leads to 



mnf} 

2h^N 

= M^(T,So), (3.30) 

where we have introduced the layer magnetization 

M,(T, Bo) = Yl (^"7 + <7> • (3-31) 

n 

An analogous result is found for the Periodic Anderson Model where one also 
has to take into account two spin systems, associated with the conduction and 
the /-electrons, respectively. Thus, with a slightly modified definition for the 
magnetization (where, for reasons of simplicity, we choose identical prefactors 
gjfiB/hioT both spin systems), 

M.(T,5o) = E--^""" (<. + -L) ' (3-32) 



eqn. ( 3.30| ) holds also for the PAM. The same formula (|3.30|) , obviously, is true 



for the Heisenberg and the Hubbard model, the only difference being the fact that 
in these cases only one spin system contributes to the magnetization. Eqn. ( p. 301 ) 



is therefore valid for all the four models discussed in this section, as indeed one 
would expect from the very general considerations that led to our choice for the 
operators A and C. 

The RHS of the Bogoliubov inequality ( p.21| ) is proportional to the anticom- 
mutator sum 

(3-33) 

k 

which, for the s-f model, is given by 



kmn 



(3.34) 



Summing the exponential over all k gives the delta function N6mn- The mixed 
commutators involving both spin and creation/annihilation operators in this case 
do not vanish, since we have the anticommutator rather than the commutator. 
Thus, we arrive at the expression 



+ {^n'yl^n^^C'^^lCn-yi + ^n^yl'^n^l'^n-fl^n-y^)) ■ (3.35) 



Since we are interested in the quantity [A, as an upper bound in the 

Bogoliubov inequality, it suffices to find upper bounds for the individual expecta- 
tion values on the RHS. For the Heisenberg term we find 



E ( fc' = 2 E + (s^n^r) < 2 e {^n,} = ^^'s{s + m 

n n n 

(3.36) 

For the Hubbard contribution we have 

n n 

n 

< 2N. (3.37) 
For the mixed terms appearing in the s-f model, an upper bound is given by 



E {4^.Cn^-aS-;) < 2 (4 + 2S{S + 1)) h^N (3.38) 

as will be shown in the following mathematical digression (see also the derivation 
of ([B.6| ) in appendix B). 



Mathematical digression: Constructing an upper bound for {^^n-yc^n-y-aSn^) 

The problem in evaluating this expectation value arises from the fact that we are 
dealing with seemingly unrelated operators. It is, therefore, convenient to use the 
identity 

„ C-cr . 





-a + (^C^7- 








Cfi'ya ^wy ) 




-a + iS^j) (^C^7- 


-a^n-ya '^^wy ) 






-cr^wya + ^'^'^7 ) } 









(3.39) 



where j labels the individual terms and is a phase factor (0(1) = +1, 0(2) = —1, 
0(3) = i, 0(4) = —i). The original expectation value has thus been decomposed 
into the sum of expectation values of pairs of adjunct operators Sjo-, -Bjo-- Using 
the spectral theorem, we find 



4 oo 
■' — oo 



(3.40) 



where 'S'^+o. {E) — ■^(\B'^^,Bj„^\{E) is the spectral density in its energy 

representation. The LHS is the expectation value of the sum of two adjunct 
operators and is therefore real. We may now make use of the fact that for pairs of 
adjunct operators, the spectral density is positive definite, which, together with 
the triangle inequality and |0(j)| = 1 gives an upper bound for the RHS: 

^ oo ^ oo 

(3.41) 



— OO — OO 



The sum on the RHS now consists of the 0-th spectral moments associated with the 
operator pairs Bj'^, Bj^,. Each of the spectral moments is given by the expectation 



value of the anticommutator 



oo 



\ j "^^^sJE^ - = ■ (3.42) 



The anticommutators can be easily evaluated: 



-Bi -Bio- \ L '^7""' " n7 ' "n^ya""'!^'-' ' n7J-(-^ 

nn^_„ + rin-y-a) ^ +2S{S + l)h 

-\-2h (^S^^C„^_^Cn-ya + »S'„^ C^^^Cn-y-a') 
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+ 23(5 + l)h'' + 2h (S:^c+^_ 






B2„B2a 
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Ah^ + 23(3 + l)h^ - 2h{3^^c+^_ 


_|_ Q-'^r^ r 
a-^n-ya- ~ ^n-y '^n-ya^n-y- 
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B^^Bs^ 
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Ah^ + 23(3 + l)h^ + 2M (5:^c+^ 


r — r 

-a^n-ya ^n'ycr^'n-y 
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Ah^ + 23{3 + l)h^ - 2hi (5^^c+^ 


r — r 

-cr^wyu ^wycr^n-y 




and hence 











J] (c+,,c„,_<,5-;> < 2 (4 + 2^(^ + 1)) (3.43) 



which proves the upper bound ( p.38 



An analogous treatment of the PAM anticommutator sum leads to a similar 
result, so that one can tabulate the results for (^[A, ^'^] + ) get... 

• ...for the Heisenberg model: 

^(^[A,A+]_^^ <23{3+l)h^N^ (3.44) 

k 

• ...for the Hubbard model: 

^{[A,A+]\<2h^N^ (3.45) 



...for the s-f/Kondo-lattice model: 

J2 ( [A A+] ^ < {4S{S + 1) + 10) fi'N^ 



(3.46) 



.and for the Periodic Anderson Model: 



(3.47) 



3.2.3. The double commutator {[[C, H]_,C+]_) 

The double commutator (^[[C, H]_,C^]_) still remains to be calculated. For rea- 
sons of simplicity, we shall, in this section, describe the calculation for all three 
models step by step. The Hamiltonians all decompose into 

H = Hq + Hf,, 

where Hb denotes the contribution to the Hamiltonian due to the external mag- 
netic field. 



The Heisenberg case 

The double commutator ([[C, H]_, C"'"]_) in this case is 

[[C,H]_,C^]_) = Y.(\[s;^{k),Ho + H,]_,S:{-k)\ \ 

\ L J — / 



(3.48) 
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The real-space commutator on the RHS is easily evaluated to give 

\C+ ff] C- _ _plpP'(AQz QZ I C- C+ _|_ C+ C-\ 

L m7' J _ ' pe pm \ pe m-y ' pem-y ' rwy pe J 



+2bh^e 



ca 

2-iK-Rr, 



(3.49) 



Inserting this into the full expression (|3.48|), we arrive at 



•ye mp 

' {'^Sps^m-/ + ^mySpe)) (3.50) 



£ m 



To this we add the double commutator {[[C, H]^, C^]_) i~k) which, as discussed 
above, is a positive real number. Replacing the spin operator expectation values 
by the upper bound 2h^S{S + 1), we find, after some minor algebra. 



4 
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(3.51) 

where we have already used the fact that 

1 - cos(fc ■ (i?„ - R^)) < ^^^^-^-^p)\ (3.52) 

With the above definition of the constant Q, we arrive at the final result 

{\\C,H]^,C+]_^ < ANdff (^\BoM{T, Bo)\ + 2h'^S{S + l)Qpy (3.53) 

The Hubbard case 

Again, we need to calculate the full double commutator (^[[C, H]_,C^]_y With 
the Hubbard Hamiltonian given in eqn. ( p.6|) , we have 



a 

(3.54) 



which still needs to be commuted with (~^)- Replacing the expectation 

values by their modulus, one gets the relevant inequality 



7£ 



a+ik),Ho + H, 



-k) 



ilae 

\{cUc,^l)\) +2hh'Y.{<^^)) ■ (3.55) 



For the same reasons as above, we may now add ([[C, H]_, (— /c) to get an 

upper bound for the LHS. The expectation values | (c^^-fQet) | + Kcztj,Cjaj.)| can 
be bounded from above by a procedure analogous to the one carried out in eqns. 
( P.39| )-( ^I4^ ), leading to |(c+^Qe|)| + < 4 (see appendix B, eqn. (p^), 

for detailed calculation). In total we may thus write, using the notation above. 



[[C, C+] < ANdh^ {\BoM{T, Bo)\ + 2gP 



(3.56) 



The s-f case 

The s-f double commutator 



[iC HU C1 _) = E ( f K(^) + ^0 + H,]_, s;{-k) + a;{-k) 



■ye 



(3.57) 



may be computed by considering first the field-independent contribution to the 
Hamiltonian 



and the respective commutators 



(3.58) 



(3.59) 



[^i, ^o] _ = Jh - 25^^^+^) . (3.60) 

We then find (neglecting for tlic moment tlie external contribution Hb) 

kmn f^-yS 



kn p'yS 

■ {W^4siCk-yi + Wt4p^Cr,S^) h\ (3.61) 



Following the usual procedure of adding ([[C, C+]_) {—k) we arrive at the 
upper bound 

{[[C,H,UC^]_) < 2fi^5^|r„t|(l-cos(^-(^„-4)))- 

•(l(4mCfc7i>| + |(c^/3TCfc7T>|) 

= 2Ndqh^k^, (3.62) 

where |(c^^j^Cfc^j.)| + | (c^/3|Cfe^T) | — already been used. 

We still need to add the double commutator ([[C, Hb]^, C+]_) with the exter- 
nal part Hb of the Hamiltonian to this preliminary result. Here, too, we need to 
take into account the presence of two distinct spin systems. Thus, Hb is given by 

= -^0^ T.(9JSL + 2aL)e-^^-^% (3.63) 

ia 

where the usual definitions of the spin operators apply. As in the previous cases, 
one gets 



[[C, Hb]-, C+] _j = -2h^NdBoM{T, Bo), (3.64) 
where, however, the magnetization is now defined as 



Finally, we arrive at 



[[C,H]_,C+]_^ < 2Nd!f (^\BoM{T,Bo)\+qP^ . (3.66) 
The PAM case 

For completeness, we also present the results for the Periodic Anderson Model. 
The detailed calculations do not differ much, in spirit, from the ones carried 
out above for the Hubbard and the s-f model, except for the difference in the 
(anti) commutation relations ( |3.23|) as compared to the s-f model. The external 



part Hb of the Hamiltonian again leads to a contribution to the double commutator 
that is proportional to the total number of lattice sites N ■ d, the external field 
Bo and the magnetization M{T,Bo), and its expectation value is bounded from 
above by 

[[C,i/,]^,C+] ) < ANdh'' \BoM{T,Bo)\ . (3.67) 

What is noteworthy about the double commutator with the unperturbed PAM 
Hamiltonian Hq is that only the hopping part Hg leads to a non-vanishing con- 
tribution, while Hf, Hff, and Hhyb do not even "survive" the first of the two 
commut at ions ,0 

[C, Hf]_ = [C, Hff]_ = [C, H,yb]. = 0. (3.68) 
Thus, the double commutator is essentially of the Hubbard form 

<ANdff(^\BoM{T,Bo)\+2qpy (3.69) 

The fact that the bound on the double commutator is independent of the hy- 
bridization and other parameters in the Hamiltonian, except for the hopping, 
confirms similar results obtained for the purely two-dimensional case fST 



^The fact that in the PAM the anticommutation relations ( 3.23 ) hold (instead of the commu- 
tation relations for the s-f model) is irrelevant for the commutator with Hs, Hf, Hff, as in these 
cases only pairs of c- or /-operators occur. For the hybridization part Hhyb it is straightforward 
to show that in this case, too, already the first of the two commutators vanishes. 



3.2.4. Proving the absence of spontaneous magnetization 

Within aU the models discussed (Heisenberg, Hubbard, s-f/Kondo-lattice, and 
Periodic Anderson Model), we have found for the double commutator 

( [[C, H]_, C+] _) < eoNd {\BoM{T, Bo) \ + ^P) , (3.70) 

where are constants depending, at most, on fixed parameters of the respective 
many-body models. We also know that 

{[C,A]_) = ^2NM,{T,Bo) (3.71) 

and 

J2([AA^]\<2^sN'. (3.72) 



We can now give a generally applicable discussion of the layer magnetization. 
With the Bogoliubov inequality ( p.21| ). 



fc ' ' k 



we find 



^ e^N'^M^(T,Bo) 

eoNd(\BoM{T,Bo)\+^ik^^ 



The LHS of the inequality can be replaced by an integral, using the formula 



j:~ld% (3.75) 



k 



where is the area in two-dimensional fc-space associated with one quantum 
state. Restricting the support of the integral to a finite-volume sphere inscribed 
into the first Brillouin zone only strengthens the inequality, so 



1 L\^,^ 



ko 



where ko is the cutoff corresponding to the sphere in /c-space. In the thermody- 
namic hmit, L and are taken to infinity (L, N — > oo) in such a way that the 
two-dimensional specific volume v^^^ = ^ approaches a finite constant. Evaluat- 
ing the integral and performing some minor algebra, we have 



again is an unsignificant constant). 

From this formula, it is clear that, for any finite temperature (/5 < oo) and 
finite thickness {d < oo), the logarithm in the denominator will diverge in the 
limit ^0 0, thus forcing the layer magnetization My to vanish. One should 



both ferromagnetic and antiferromagnetic ordering. Thus, we have proved that 
the Mermin- Wagner theorem which was originally shown to exclude a magnetic 
phase transition at finite temperature for one and two-dimensional systems can be 
extended to Heisenberg, Hubbard, s-f (Kondo- lattice), and PAM films of any finite 
thickness. While this was to be expected from very general considerations, we have 
been able to confirm this conjecture through a detailed microscopic calculation 
that takes into account specific Hamiltonians and distinguishes between individual 
layers. 

Our calculation suggests that the parameter d, i.e. the number of layers in 
a film, plays a similar role as the inverse temperature /3. If either one diverges, 
a phase transition cannot be ruled out. It is, of course, not a priori clear if the 
behaviour of an upper bound on a physical quantity has any physical significance 
itself; however, it is reasonable to take this result as indicating that in order to 
describe (anti)ferromagnetism, one has to take the thermodynamic limit seriously, 
in the sense that N ^ oo does not suffice, but d — > oo is required as well. 

3.3. Problem of superconducting long-range order 

The disovery of high-temperature superconductivity with its surprising character- 
istics has opened up a new field of research within many-body theory. The details 
of the origin of superconducting long-range order in the high-Tc superconductors 
are still controversially debated, but it seems certain that the transition to the 
superconducting state is due to the strong electron-electron correlations. Thus, 



M^{T,Bo)<i 



(3d 



(3.77) 





the same set of quantum many-body models used for describing spontaneous mag- 
netic order in itinerant-electron systems is also appropriate for the superconduct- 
ing phase transition. In order to demonstrate that our discussion of the absence 
of a magnetic phase transition in film systems can be extended to cover the super- 
conducting transition as well, we shall carry out the calculations for Hubbard films 
of finite thickness which are believed to reflect the main properties of electrons 
moving in the set of coupled CuO planes which is characteristic of many of the 
high-temperature superconductors under investigation. 

3.3.1. s-wave Cooper pairing and generalized rj pairing 

In superconductors the essential new feature is the formation of Cooper pairs, i.e. 
of pairs of conduction electrons with opposite spin and opposite wave vector 

which travel freely through the crystal lattice. Before one attempts to apply 
the Bogoliubov inequality to investigate the possibility of a phase transition to 
the superconducting state, one must decide on a quantity, or order parameter, 
to characterize the phenomenon. For (anti)ferromagnets this decision was fairly 
obvious, since the macroscopic quantity measured experimentally, i.e. the (stag- 
gered) magnetization, is straightforwardly given, on a microscopic level, by the 
sum over the 2;-component of the spins. Instead of characterizing superconductors 
by their vanishing resistivity on the macroscopic level, we shall focus on the un- 
derlying microscopic process of pairing between electrons, which is characterized 
by a breakdown of U(l) symmetry of the conduction electrons. 

Before giving the Hamiltonian for Hubbard films with an external U(l) sym- 
metry breaking term, we briefly summarize the local on-site pairing operators 77^ 
which describe the creation or annihilation of a Cooper pair at lattice site {i, a) 
where i denotes the position in the 2D Bravais lattice and a as usual is the layer 
index: 

Vta = c+T^-i (3-78) 

Via = {Vta}^ 

= CialCia]. (3.79) 

We also introduce an operator r/^'^deflned as 

Via = \ i^ia^ + ^ial " 2) . (3.80) 



Going over to fc-space by Fourier transformation, one can verify the commutation 
relations 



(3.81) 
(3.82) 



With ( p.78|) and (|3.79| ), a measure for the breakdown of U(l) symmetry due 
to pairing is given by the expectation value^ 



77+(-i?) 

N 

N 



(3.83) 



For film systems (and since we are concerned here with local on-site pairing), we 
can, as in the magnetic case, write the total order parameter ^{K) as the sum 
over intra- layer quantities daiK) : 



a i 



(3.84) 



The ordering wave vector K allows one to distinguish between different types 
of pairing: for K = one would have s-wave pairing, for K (generalized) 



//-pairing ||39| , 

The Hubbard Hamiltonian with a U(l) symmetry breaking contribution of 
order A is then given by 



H 



ija/3a 



lOLU'^lOL — CJ 



-iK-Ri 



-iK-Ri 



(3.85) 



*We follow the notation used in similar proofs for the purely one- and two-dimensional cases 

([io|,ii). 



where the second hne is the symmetry-breaking part 
In order to make use of the Bogohubov inequahty 



E 



(3.86) 



to derive an upper bound for '$a{K), the operators C and A must be chosen 
properly. With 



7 



(3.87) 
(3.88) 



one finds that indeed the Hamiltonian-independent commutator [C,A]_ repro- 
duces the layer-dependent quantity da{K) : 



{[C,A]_) = ($^[r/^(fc),r/^(-fc-i?)^ 



(3.89) 



The anticommutator sum + ) t)e shown to be bounded from 

above: 



k V 



(3.90) 



The double commutator, after some algebra, can be reduced to the familiar 
structure quadratic in k plus a quantity proportional to the order parameter: For 
the contribution due to the unperturbed Hamiltonian Ho (i.e. the first line of 
( ^.85| ) only), the double commutator is 



< NdQk'^ 



(3.91) 



7? 



where, as usual, J^im-ye '^im {^i ~ = < oo is assumed. The double 

commutator with the U(l) symmetry breaking term Hgc (second line in (|3.85|) ) 
is bounded by a quantity proportional to the total number of lattice sites, the 
total order parameter ^{K), and the parameter A with which the U(l) symmetry 
breaking can be turned on (A 7^ 0) or off (A ^ 0) : 



7£ 



< 2XNd 



(3.92) 



The absence of s-wave (or generalized rj) pairing in Hubbard films of finite 
thickness now follows straightforwardly from the Bogoliubov inequality ( p. 861) in 
the thermodynamic limit. Formula ( p.86[ ) with ( |3.89| ), ( p.90[ ) and the sum of ( p.Ql] ) 
and (|3.92|) , is equivalent to 



E 



Ar2 




2 


NdQk^ + 2\Nd 





< AN 



(3.93) 



so that in the thermodynamic limit, with the familiar substitution 

J2- ^ {V^''^/{27tY) J d^k..., 



we get for films consisting of d two-dimensional layers stacked on top of each other, 
the result 



,(2) 



27r 



< 



Pd 



In 1 + 



2\\d(K)\ 



(3.94) 



which in the limit A — > proves conclusively that no pairing transition of the 
proposed kind and, thus, no corresponding superconducting phase transition can 
occur in the Hubbard model, provided (3 and d are finite. 



3.4. Two remarks on special cases of low-dimensional order 



3.4.1. Planar magnetic order 



Ferromagnetic and antiferromagnetic ordering along the z-axis are not tlie only 
types of ordering conceivable. Indeed, for the Hubbard model so-called spiral 
states have been discussed, which are conjectured to play an important role in 



the thermodynamics of the two-dimensional model at low temperatures In 
particular it has been suggested that in the zero-temperature case at half filling the 
2D Hubbard model with U = oo may display an instability towards the formation 



of spiral phases [E2 



In this context, it is of interest to consider order parameters which keep the spin 
vectors confined to a plane. It has been shown that any such order parameter 
must vanish at finite temperatures in ci < 2 dimensions. By analogy with previous 
expressions for the magnetization, a planar magnetization MpianiT, Bq) is defined 
as 

MpianiT, B^) = ^Y. + /^^O) ' (3.95) 

i 

where a^, (5i are arbitrary real numbers except for the requirement that the sum 
of their squares be convergent: 

^ E H + (^) ^Q<N. (3.96) 

i 

With the Bogoliubov inequality ( |2.21| ) and by choosing A, C so that 

A{k) = E e-^'-"^^ (/?,aj - «,a|) (3.97) 
j 

and 

B{k) = Y,e~''^'^'(T], (3.98) 

j 

it is easy to verify that for planar magnetic order, too, the order parameter, 
Mpian{T, Bq), can be bounded from above (^,^i = const), 



so that in the hmit of vanishing external field at any finite temperature /3 < oo, 
the planar magnetization vanishes: 



In this remark, we have sketched an extension of the Mermin- Wagner theorem 
to a case of magnetic order where ordering within a plane perpendicular to a fixed 
direction is considered. The details of the calculation are given in we shall 
return to this case briefly in a later section of this thesis where it serves as a simple 
example for the application of a more general correlation inequalities approach. 

3.4.2. Fractal lattices with spectral dimension d <2 

We shall briefly outline another interesting special case of the Mermin- Wagner 
theorem, namely its derivation for fractal lattices whose spectral dimension is 
d <2. In the previously discussed cases, extensive use was made of translational 
invariance of the underlying 2D Bravais lattice in order to exclude a magnetic (or 
pairing) phase transition in film systems. Spatial dimensionality, in these cases, 
comes into play in the thermodynamic limit, when the sum over wave vectors 
is replaced by a two-dimensional integral in fc-space. Thus, the translationally 
invariant lattice is naturally embedded into a Euclidean space of integer dimen- 
sionality D which determines the continuum limit of the statistical system under 
consideration. The question, however, remains how this embedding into Euclidean 
space can determine the system's behaviour when no explicit reference to its di- 
mensionality is contained in the lattice Hamiltonian. A detailed analysis of the 
mathematical physics of this problem would be far beyond the scope of this thesis 
(see [|3] for an elaborate discussion). 

It turns out that for irregular lattices, such as fractals, one must adopt a 
more sophisticated notion of "dimension" which is allowed to take non-integer 
values ("spectral dimension"). By use of a mathematical theorem on divergencies 
in Gaussian field theories (a proof of which, along with an extensive discussion, 
is given in ^^) it has been shown that also for the Heisenberg model on fractal 
lattices with (spectral) dimension d < 2 and nearest-neighbour interactions, the 
spontaneous magnetization must vanish at any finite temperatures [^. 

The extension to other situations or models, as well as the physical meaning 
of this result, are far from obvious and we must, at this point, refer the reader to 
refs. [^], for a more detailed discussion. 



4. Green's Functions and Correlation 
Inequalities Approach 



4.1. Linear response quantities and their algebraic proper- 



4.1.1. General formalism of linear response theory 

In examining the reaction of a physical system to an external, or in fact any, 
perturbation, one will generally be faced with a Hamiltonian H = Hq -\- H^xt, 
where Hq includes the kinetic (or, more generally, the one-particle) term and the 
internal interaction due to, say, a pair-potential. The external contribution Hg^t is 
due to an external scalar field Ft which couples to an observable B of the system 
and which may be explicitly time- dependent {He^t = FtB). This perturbation 
will lead to a change in the system's behaviour affecting other observables, such 
as A, as well.Q The dynamics of the process are given by the time-dependent 
Schrodinger equation 



where in the second step we have switched to the interaction picture defined by 



ties 




H\<^) 



(4.1) 



ih-\^n{t)) = Hgtit)\^D{t)) 



(4.2) 



with 



^iHot/hjj^^^^-iHot/h 



(4.3) 



^This section follows the derivation in 



i 



41 



(the substitution with Hq — Hq — jiN in the last step is possible if [Hq, 7V]_ = 0). 

Evaluating the change in the expectation value of the observable A requires 
calculating the trace of A with the grand canonical density operator at any given 
time t : 

= tr [ptA) . (4.4) 

The time- dependent density operator within the Schrodinger picture satisfies the 
von Neumann equation of motion 

ihpt = [Ho, pt]- + [-ffext, Pt]-- (4.5) 
In the interaction picture, pt transforms as 

pf{t) = exp i^^not^ Pt exp (-;^^o^) (4.6) 

so that, taking the free initial condition 

1- _ exp(-/37io) . , „^ 

into account, one has 



t 



P?{t) =P^-lj dt' [HgAt'lP?{t')] _ (4.8) 

— oo 

which can be solved by iteration (the somewhat redundant prime in " ext'" in- 
dicates that the explicit time dependence of i^ext is also to be integrated over). 
Averages with the density operator &X, t — — oo, i.e. with the field-independent 
density operator po, will be denoted by an index "0": 

i^^ = tr (poA) . (4.9) 

To first order, the density matrix thus is 

t 

pf (t) ^Po-^J dt' [HgAt'):Po] _ (4.10) 



— oo 



or, in the Schrodinger picture (using ( [4.3|) ), 



Pt ~ Po - / dt' exp 



[^e^t'(i^'),Po]_exp 



iHot 



(4.11) 



With this expression, it is now possible to exphcitly evaluate the change in A, 



A 

i 

"h 



A 



dt'FH 



A^{t),B''{t') 



(4.12) 



The integrand in ( [4.12| ) depends only on the commutator of the observables in- 
volved and on the "internal" dynamics of the interacting many-body system at 
zero external field as described by Hq, which enters the calculation via the averag- 
ing procedure ([4.9|) . To emphasize this fact, one introduces the retarded Green's 
function 



ret 



Ait)- Bit') 

-te{t-t') ( A{t),B{t') 



(4.13) 



For e = +, the commutator on the RHS is the same as in ( [4.12| ), because at 
zero external field the Heisenberg picture coincides with the interaction picture. 
Dividing the expectation value of the (anti-) commutator in ( [4.13| ) by 27r, we define 
the spectral density 



SABit,t') = ^( Ait), Bit') 



(4.14) 



which is of central importance and which will be discussed later (e.g. in ( [4.38| ) in 
its spectral representation in Fourier space). 

Thus, after Fourier transformation (which requires regularization by use of an 
infinitesimal damping factor with ir] = iO'^) 



AA, 



dEFiE)G%iE) exp [-j{E + t7])t 



(4.15) 



which is known as the Kubo formula. It can be interpreted as describing the 
response of a system to a perturbation in B in terms of the change in the observable 
A. 



Equation of motion for the Green's function 

From elementary quantum mechanics it is known that in the Heisenberg picture 
a (time-dependent) operator At{t) evolves according to the equation of motion 

ihj^Ait) = [A,n]_ {t) + ih^. (4.16) 

Calculating the equation of motion for the Green's function GABit,t'), which 
consists of combinations of operators, leads to 

th^GTsM = h6{t-t') <[A5]_,> + {{[A,H]_ {thB{t'))y'' (4.17) 
or, in terms of energy, 

E {{A; B))f = h {[A, BU) + (([A n]_ ; B))';' . (4.18) 

By iteration, a hierarchy of higher Green's functions is generated on the RHS 
which, if infinite, can be decoupled artificially (by an appropriate approximation). 
A procedure for decoupling the higher-order Green's function in the Heisenberg 
model will be presented in the next chapter. 



4.1.2. Transverse spin susceptibility in the Hubbard model 

For a system with an external magnetic field, the response of the magnetization 
is the quantity one is interested in and which one can hope to express in terms of 

— * 

Green's functions. The external contribution H^j-t = —rh ■ Bt to the Hamiltonian 
(m: total magnetic moment; Bf. time-dependent external magnetic field) will 
lead to a change in the magnetization M = y (m) . Thus, by virtue of the Kubo 
formula, the response of the (3-th component of the magnetization to Bt is given 

by 

oo 

= ^ E / (("^' "^"(^')» • (4.19) 



This formula suggest defining a magnetic susceptibility by 

xr(t,t') = ((^fW;-- (O)) (4.20) 

which generahzes in a straightforward way the more phenomenological interpreta- 
tion of the susceptibihty as a "proportionahty coefficient" between external field 
and magnetization, yielding 



AM, 



13 ^ 



-E J dt' (j24''it,t')B?] . (4.21) 



For (/?, a) = (+, — ) eqn. (^4.20|) describes how the system reacts to spin fiips. 

With (|1.12|) and ( [L.13D the transverse spin susceptibility x'^~{E) of the Hubbard 
model is given by 

^ ~Vh (4.22) 

or equivalently in /c-space 

= (4-23) 

kp 



where 

/i 

and we have introduced as a shorthand 



7 := f toB)' (4.24) 



4.1.3. Dynamical structure factor 

The dynamical structure factor C is anotlier important physical quantity which, 
as we shall see, provides the same insight into a system's behaviour as the suscep- 
tibility [^] , . We shall derive it from the transverse two-particle spin correlation 
function 

Ctit-t') = {^tm(t')) (4-26) 

by Fourier transformation: 



oo 



■iq-(R,-Rj) 



■' — oo 
oo 



1 / d{t-t;){o\q,t)a-[-q,t;))e^^'-'^l\ (4.27) 



4.1.4. Connection between susceptibility and dynamical structure fac- 
tor 

Along with ( |4.27| ) we also know, by definition, that 

- R,- t-t') = -ze{t - t') {{antyjit')) - {crj{t')ant))) (4.28) 
for the retarded Green's function and 

oo 
-oo 

- (a7(0)(T+(r))} e'^-(^»-^.) (4.29) 



Ato(5'jAti?) 



NVh 



a-{k,0)a+{-k,r)jj. (4.30) 
By use of the integral representation for the Heaviside step function, 

oo 

1 r iE'r/h 

e(r) = / dE'^ , (4.31) 



the susceptibility can indeed be expressed completely in terms of the dynamical 
structure factor: 



oo oo 



NVh 2n J ^' J ^" {E"-E)-ir] 



dr / dE' 



— oo — oo 
oo 



l^o{gjl^Bf\ }_ f ^j^,c+-{k,E")-c-+{k,-E") 



Vh J 2tt J {E- E") + ir] 



-oo 
oo 



.„toB)'^ 1 [^^.(i-^^nc^-ik,n^ ^^^^^^ 



Vh J 27: J {E- E") + iT] 

— oo 

where use has been made of the quite general relation 

Cab{-E) = e-^'^CBAiE) (4.33) 

which can be verified by inspecting the spectral representation of the dynamical 
structure factor, 

oo 

C^-{q,E) = ^ j dte'^'/^{a^q,t)a-{-q,0)) 

— oo 
oo 

1 r ^ p-l3En 

n ^ 

— oo 

oo 

= Y.^;^ / dter^^+'^^--^^^^' {n\a+{q,0)\m) {m\a-{-q,0)\n) 



= 27rn^^^|(n|a+(g,0)|m)p5(E+(E„-Ej). (4.34) 

nm 

From this expression one also concludes C^~{q, E) > 0, which is true for any dy- 
namical structure factor C^+a(^5 E) corresponding to a pair of adjunct operators 
A, A'^. Thus, the dynamical structure factor is positive semi- definite. The equa- 
tion for the susceptibility ( ^4.32| ) can easily be generalized for arbitrary operators 



A,B 



XAB{E;ir]) 



1 

2^ 



dE' 



-oo 

oo 



,,Ca b{E') — Cba{—E') 
E-E' + irj 

,{l -e-^^')CnA{E') 
E-E' + ir] 



(4.35) 



The generalized susceptibility still depends on the energy E and, on a more formal 
level, the regularization if]. Often, one is interested in the static limit, where 
E ^ 0. The static susceptibility is then defined as 



XAB = ^ lim XAB{E;iri) 



and thus, for /c-dependent operators A{k), B{k'), 



XAB{Kk')^^\ lim — 
E,'n-*o Zn 



dE' 



, (1 -e-P^')CBA{k',k-E ') 
E-E' + ir] 



(4.36) 



(4.37) 



The latter formula will later serve as the starting point for the discussion of 
certain algebraic properties of the static susceptibility that will allow us to regard 
Xab as defining a scalar product between A and B in operator space. 

At this point, however, it appears useful to connect the static susceptibility 
with yet another ubiquitous quantity, i.e. the spectral density, which in its spectral 
representation is {e — ±1) 



h 



Sil{k,k';E')^-Y,{n\A{k,0) 



m) (m 



B{k',0) 



n > e 



-PEr, 



[e^^' - e) S{E' - {E„ 
(4.38) 



Comparing this with the spectral representation of the dynamical structure factor 
CBA{k\k-E') 



n 



(4.39) 



oo 



-I3E. 



NE 



n) (n 



m ) e 



J{E'-{E,n-E^))t/h 



m 



A{k,0) 



n)(n 



B{k',0) 

B{k',0) m)6{E' -{Em-E,,)) 



and reminding oneself of ([4.35|) , it is easy to see that in the static limit 



'^i )iT^^XAB{k,k';E,ir]) ] = XAB{k,k') 

\E,ri^O 



•/^A sikj k'; E') 
E' 



V / dE' 



(4.40) 



Algebraic properties of the static susceptibility 

Earlier it was mentioned that the algebraic properties of the static susceptibility, 
notably the fact that it defines a scalar product in the space of operators A, B 
lead to interesting exact relations between different physical quantities.^ Three 
properties need to be verified in order to prove that 



(4.41) 



A{k)-B{k')):=XAB{Kk') 

is a valid scalar product: 

1- XAsik, k') is linear in A{k) and B{k') 
2. {xAB{k,k')X =XB^A+{k\k) 



3. XAA+iJ^) = XAA+ik, k) = V J dE' 



„S^^+(k,k;E') 



> 0. 



^This section partly follows ref. 



The linearity is trivial, as can be seen directly from the definition of Xab- The 
last two properties can be demonstrated by use of convenient representations in 
terms of either the spectral density or the dynamical structure factor. For 2.) we 
start from eqn. ( [4.35| ): 



oo 



2tt J {E-E' + ir])* 



oo 

^-P I dE 



271 J E-E' -ir] 

— oo 



oo 

—V I dE' 



2n J E-E' -ir] 

— oo 

= XB+A+{k',k;E,-ir]) (4.42) 

so that in the static limit, where £H (limg ,,^0 •••) is to be calculated, one indeed 
has 

[xAB{k,k')y = XB+A+{k',k). (4.43) 

For property 3.) it is easiest to start directly from the integral representation 
of the static susceptibility (^4.40|) : 

XAA+ik,k) 

= XAA+ik) 



oo 



E 



-oo 
oo 



= V J rfi?l|5^|(n|A+|m>|^-/5^"<5(i?-(i?„-i?J)(e^^-l) >0 

mn 

— oo 

by use of the spectral representation for the spectral density in the last step. Thus, 
the static susceptibility indeed defines a positive semi-definite scalar product in 
the operator space. 



4.2. Correlation inequalities approach 



Having shown in the previous section that the static susceptibihty defines a scalar 
product between A and B, we shall discuss, on quite general grounds, an inequality 
for the static susceptibility. The inequality, in fact, will be a special case of the 
Schwarz inequality which holds for any scalar product and, thus, in the case of 
the susceptibility reads 



XABik.k') <XAA+{k)xB+B{k') 



(4.45) 



For the specific example, we require A^{t) to be of the form A^{t) := 

The spectral density for the "original" operators A and B can then be expressed 

in terms of B and the newly introduced operator Q : 



SABik,k';E) 



d{t-t')er^''^'-''^SAB{k,k';t,t') 



oo 



E(t~t') 



dQ^jt - t') 
d(t - f) 



(4.46) 



The next step is to integrate out the RHS by partial integration. For any physical 
problem, one can multiply the integrand with a factor e~"l*~* a > 0, and take 
the limit lim^^o after the integration has been carried out. Thus, contributions 
from the infinitely distant past or future are neglected. In this case, one arrives 

3See e.g. ref. Isotl. 



at 



SAB{k,k';E) 



1 

2^ 



-oo 
oo 



l-E J d{t-t')e^^^^~^'^{[Q^,{t-t'),B^„{0)]_ 



27T 

= ESQB{k,k';E). 
With an analogous calculation it can be shown that 

SAAAk;E) = E^SQQ+{k;E). 



(4.47) 



(4.48) 



With these formulas one can now make use of Schwarz's inequality for the 
susceptibilities ( [4.45|) . Rearranging the latter as 



XABik,k') 
XB+B{k') > 

XAA+{k) 

and deriving from ( |4.40|) the expressions 

XAB{k,k'-t,t) = ([Qj^{t),Bp{t)] 



(4.49) 



(4.50) 



XAA+ik) 



one has 



/ dESsMk';E) ^ 

E 



dEESQQ+ik,E) 



(4.51) 



/ dEESQQ+{k,E) 



(4.52) 



Since the LHS of this inequahty is essentially identical with XB+B{k'), eqn. 
( [4.52| ) defines a relation between the response function for the observable B and 
commutators which can, in principle, be calculated directly from our knowledge 
of Q, B, and H. 



4.3. Finite-temperature case and Bogoliubov inequality 



So far in this chapter, quite general relationships between different physical quan- 
tities have been stated. In order to show that the formalism developed up to this 
point is also capable of reproducing the substantially stronger results obtained in 
Ch. Ill by use of Bogoliubov's inequahty, it is useful to verify that the latter can 
indeed be obtained within the framework of this chapter. 

Pitaevskii and Stringari have suggested to define a scalar product simply 
through the anticommutator: 



{A;B) ■.= {[A\B]^) 
which gives the Schwarz inequality 

{[A\AU){[B\BU)>\{[A\Bl 



(4.53) 



(4.54) 



It is possible |^ to define auxiliary operators, denoted by a tilde, in a such a way 
that 



n 







n 



from which it follows that 



n 







n 



(n|qo) 

-(0|C|n) 

-(n|C+|0) 
(0|C+|n). 



(4.55) 



(4.56) 



If we assume, for the moment, that the expectation values in (|4.54| ) are evaluated 
in the ground state (T = 0) and instead of B the newly defined operator B is 
used, one deduces from 



(0|[A+,5]_|0> 
(0|[5+,5] + |0> 



(4.57) 



that the Schwarz inequahty holds not only for the anticommutator on the RHS 
but also for the commutator, so that (writing down (|4.54| ) for the operators A, B, 
and making use of ( 4.57|) to express B in terms of B) one finds 



[[A\AU){[B\BU)>\{[A\B]^)\ 



(4.5^ 



The extension of this method to finite temperatures is simple. The modified 
operators required for this purpose (denoted by a hat) are given by 



n 



C 



m 



[n \C\ m) 



(4.59) 



where Wm = ^ exp {—PEm) is the statistical weight in the grand canonical en- 
semble. The expectation value can be calculated by expanding into eigenstates: 



[A\Bl 



mn 
mn 

{[A\B]_). 



e - e' 



■^^") (m U+ 1 n) ^"^ ^" {n\B\m) 
"^^") (m n) {n \B\ m) 



(4.60) 



If both commutators are constructed according to ( |4.59| ) , one finds 
1 



[A\Bl 



-I3E„ 



^—l3En ^—I3E„ 

fD—pEn _j_ Q—fiErr 



\ A'^\ {n \B\ m) 



mn ^ ' 

■ (m\A'^\n) {n\B\m) . (4.61) 

Once again making use of the spectral representation of the spectral density ( |4.38| ), 
one sees 



[A\BY 
or, if specifically A = 13, 



dE tanh 



= - / rfEtanh ( ^ ) Sa+a{E). 



(3E 



(3E 



Sa+b{E) 



(4.62) 



(4.63) 



This compares with the anticommutator for "normal" operators {A, B instead of 
A,B), for which 

A] + > = \j dE coth (^^^ Sa^a{E). (4.64) 

Furthermore, from the definition of the spectral density it is clear that for the 
commutator 

([A+,5]_> = dESAME). (4.65) 
The Schwarz inequality ([4.54|) may then be written as 



J dESA+A{E) coth (^(^^ j dESB+B{E) tanh 



PE^ " ^ 



2 



> 



dESA+B{E) 



(4.66) 



Earlier (Ch. II) we had encountered the definition of the Bogoliubov scalar prod- 
uct between two operators: 

En^Eui TTT" TT/" 

^iA;B) = (^l^+h) {m\B\n) 

„ En Eqji 

mn 

^ « En^Em 

= - dE (n|A+|m) (m|5|n) (e~'^^™ 

mn 

Hence, the Schwarz inequality for the Bogoliubov scalar product can be expressed 
in terms of the spectral density alone. 



From the series expansion of coth the simple relation 



(4.68) 



may be used to strengthen inequalities with appearing under the integral sign, 



'1 rdEpE 



(4.70) 



With the energy-weighted sum rule, which relates the first spectral moment 
to the double commutator 



dEESA+c{E) 



(4.71) 



and by choosing B := [H, C]^, we get the following result: 



. (4.72) 



The two remaining integrals in ( [4.72| ) are calculated in appendix C and shown to 
lead to the Bogoliubov inequality: 



(4.73) 



4.4. Low-lying excitations 

Having seen how certain algebraic properties impose constraints on the behaviour 
of the many-body quantities discussed in this chapter and having demonstrated 
that the Bogoliubov inequality can be reproduced within this framework, we point 
out, in this section, some additional applications of susceptibility and correlation 
inequalities. We do not intend to give a complete survey, as this would be far 
beyond the scope of this work, but rather outline some aspects of low-lying exci- 
tations of a many-body system and the Hubbard model in particular. 



4.4.1. Relationship between magnetization and transverse susceptibility 
at zero temperature 

The definition of the dynamical structure factor C^^ {q, E) as the Fourier trans- 
form of the transverse spin-spin correlation function still holds at zero tempera- 
ture; the correlation function, however, must now be evaluated with the T = 



ground state: 

oo 



-oo 
oo 



27Th 
1^ 



J2 - {E^ -Eo))\ (0 \a+{q, 0)\m)\ 



(4.74) 



Since Eq is the ground state energy, we have E^ > Eq, so the 5-function only- 
gives a contribution if £" > 0. As a corollary 



If one defines 



EC-{q,E)^ II iff >0 



C±{q,E) = ^{C^-(q,E) + C-^i-q,E)) 



[a+{q,E),a-{-q,E 



and 



C^{q) = j dEC^{q,E), 



one can make use of ( [4.5^ ) for operators (T+(g + Q), cr (— g) to get 



C^{q + Q)C^{q)> M{Q) 



where the magnetization 



M{Q) 



2h 

iV 



(4.75) 



(4.76) 



(4.77) 



(4.78) 



(4.79) 



has been introduced (the reciprocal lattice vector Q as usual accounts for possible 
antiferromagnetic order). Eqn. ( ^.75| ) tells us that C±{q,E) is nonvanishing only 



for E > 0, where it is positive, so the integral quantity C±{q) may be formally 
written as 



C^iq) = / dE^EC^E) 



C±{q,E) 



E 



(4.80) 



In this form, we can apply Holder's inequality (|D.2D which states that for real- 
valued function / and g for which and are integrable (where p and 
q are numbers that satisfy the relations + q~^ = 1 and p > 1) the following 
inequality holds: 



dxf{x)g{x) < i dx \f{x)\ 



i/p 



\9{x)\' 



(4.81) 



For p = q = 2 and with /, g defined as 

fix) = ^ ECME ) 
C±{q,E) 



E 



Holder's inequality ( [4. 811) applied to ( [4.80| ) gives 



Cx(g1 <^J dEEC^{q,E) jdE 
Earlier (Ch.III) we had calculated 



E 



(4.82) 



(4.83) 



dEEC±{q,E) 



a+iq,E),[n,a-i-q,E)]_ 



f{q)Nq^ < const ■ Nq^ 



which suggests writing 



JdEEC ME) _ .2,„-^ „2 



(4.84) 



(4.85) 



thus defining the new quantity uj{q). Inequality ( [4.80| ) can be strengthened using 
( [4.83| ) and rearranged by use of (|4.85| ) to give 



C±{q + Q) > 



M{Q) 



(4.86) 



This equation relates the order parameter M, the transverse susceptibihty 
and the structure factor (via C_l). The relation between these quantities can 
be made even more pronounced by summing both sides of ( 4.86|) over the first 
Brillouin zone. Drawing on our previous results for the Hubbard model we have 



\m{Q) 



(4.87) 



or, in the thermodynamic limit, 



(Tq 



M{Q) 



l.B.Z. 



q^{q)x^~{q)' 



(4. 



Provided that uj{q) and x^^io) remain finite in the limit g — 0, this means that in 
ID not even at zero temperature a transition to a state with nonvanishing M{Q) 
can occur, as this would lead to a logarithmic divergence of the RHS and thus 
would contradict being a constant. For ui{q), our previous results indicate that 
it is bounded from above. That will also remain finite under quite general 

conditions can be made plausible starting from the zero-temperature formulation 
of iWM) - 



{k,E') 



E' 



(4.89) 



where {e = —) indicates that we are now dealing with the anticommutator spectral 
density. Making use of the spectral representation one shows that at T = : 

X+"(?1|r=o = 27r/iP j dE'^J2^(^'-(^rn-Eo))- 

-oo 

|(0|a+(g,0)|m)|^ + |(0 |ff"(-g, 0)| m) |^) . (4.90) 



Integrating the 5-function gives a factor l/{Em — Eq) appearing in the sum 
over the complete set of states \m) . Thus, for \m) = |0), i.e. E^ = Eq, one will 
get a divergent term, so that at first sight this could be thought of as indicating. 



via the divergence near g = of the RHS of ( [4.88| ), a phase transition to a state 



of finite M. This divergence, however, is spurious and does not usually indicate 
a phase transition. The "excitations" which the divergent term in ( 4.90| ) seems 



to suggest, are just other ground states and thus should not be included in our 
description of a possibly ordered phase. Note that this argument should merely be 
taken as giving a hint as to the peculiarities of zero-temperature phase transitions. 
A detailed and mathematically rigorous discussion would require more elaborate 
limiting procedures as g — >■ and E ^ 0. (For related problems see a series of 
papers by Lange [i^l'EHO 



4.4.2. Goldstone's theorem 

Returning to the finite-temperature case, we want to conclude this chapter by 
determining heuristically, whether in the limit A; — > 0, assuming a ferromagnet- 
ically ordered phase, there exist elementary excitations with vanishing energy 



E(k) 4; as this would mean that Goldstone's theorem is satisfied. We start 

K — > 

from the inequality 

X^-{k) < ^(^[crHk)-{crHk)),cr-{-k)-{a-{-k) 

A J2 {d.c^.d.cr., + cl,cn,cl^,c^,) (4.91) 



mn 



which holds for any k. The inequality becomes only stronger when summing the 
LHS over all k with /c > /cq > (i.e. excluding only the null vector) while 

summing the RHS over all k in the first Brillouin zone. One then has 

mn 

Id N N"^ 

^ («)') = ^^^^^-^7^' (4-92) 

n 

where = | (c^Qi — c^Qj) has been employed. 

In the thermodynamic limit, the specific volume approaches a finite limit, 
V/N v^'^\ and the sum can be replaced by an integral. The lower cutoff ko 



becomes infinitesimal, while for any finite (3 the integral remains bounded from 
above by a finite constant: 

j d^kx^~{k) < const- 13. (4.93) 

|fc|>A;o^O+ 

From this we easily read off that as ^ 0, behaves as 

X^- - (4.94) 

with < e < 1 and / not necessarily isotropic in fc-space. 

With this asymptotic behaviour of x^" iii mind, we shall now discuss Gold- 
stone's theorem. As a measure for the mean energy of elementary excitations, we 



make use of the following modified definition of spectral moments (see also [49 
H, 0, il): 



j dEE^'-^SB+Bi.k.E) 

i?S(fc) := hm hm (4.95) 



(Note that this definition fulfills all the intuitive requirements for a mean energy, 
such as Elik) - (EBik)^ > and E'^ik) > 0.) 

The mean square value of the energy (i.e. the second modified spectral mo- 
ment) is bounded from above, via eqn. ( [4.52|) and the spectral moments, by 



Elik) < ^ (4.96) 

(where ^0 is a constant and the limit Bq ^ Q has already been taken). We also 
know that 

M^(T, 5o) < 2v^^^q ■ x^'ik) ■ e. (4.97) 

Multiplying both inequalities and taking the limit — 0, which allows us to make 
use of eqn. (|4.94| ), we arrive at 



k'° 

^'^^^ - ■ m\t^B,) W-^ =9\^\ k'--. (4.98) 




(Note that in this section we imphcitly assume that a phase transition to a state 



immediately concludes that gapless excitations exist, with the energy spectrum 
vanishing as ~ fc^~2 as A; 0. This appears to indicate that, indeed, Goldstone's 
theorem holds for the Hubbard model.^ 



^Note that our discussion is not meant to be a rigorous proof, as this would require further 
refinement of the argument. 



with finite magnetization has already 




5. Relevance of the Mermin- Wagner 
Theorem and its Validity in 
Approximative Theories 



The statement that "magnetic long-range order in low dimensions is destroyed at 
any finite temperature due to collective long-wavelength excitations" has become 
common wisdom amongst theoretical physicists. While the absence of long-range 
magnetic order is proved conclusively via the Bogoliubov inequality, it does not 
provide a detailed picture of how the destruction of order comes about. With 
respect to approximative treatments of many-body models such as the Hubbard 
model, it will be helpful to improve one's understanding of the microscopic pro- 
cesses that determine the existence, or absence, of a phase transition. In general, 
this will be a rather difficult task, due to the diversity of possible approxima- 
tions. In the following, however, we shall show that simple improvements of 
mean-field theory, though physically appealing, generally do not suffice to en- 
force the Mermin- Wagner theorem, whereas in the Heisenberg case the Tyablikov 
method already seems to incorporate the effect of dimensionality on long-range 
order in all the relevant limiting cases. 



5.1. Onsager reaction-field theory 

Reaction-field theory as an improvement of ordinary mean-field theory goes back 
to Onsager's treatment of a system of dipoles Onsager suggested that the 



orienting part of the mean-field acting on a given dipole must not include that part 
of the contribution from the dipoles in the vicinity which is due to their correlation 
with the given dipole under consideration. We shall not give a detailed review of 
Onsager's original work for a system of dipoles, but instead introduce the main 
concepts at first for the Ising model and then for the Heisenberg model. We do 
not claim any originality for the technical calculations which for the Ising and the 
Heisenberg model have been carried out by and respectively. 



63 



5.1.1. The Ising model 



As the purpose of this section is to introduce the main ideas, we shaU restrict 
our attention for the moment to the unmagnetized phase of the Ising model. The 
total mean field acting on a spin (in the sense of Weiss mean-field theory) is 
given by 

BY'''' = Br + Y.V,,{^^,), (5.1) 

where Bf^^ is an infinitesimal external field and —Vijfiifij is the spin-spin interac- 
tion between sites i and j, which must be summed over all j ^ ito yield the total 
mean field. The mean field has an orienting field on the magnetic moment of the 
spin i under consideration, thus giving 



)Weiss 



" (5.2) 



where Cq is a constant accounting for proper physical dimension and V^j = as 
usual. By Fourier transformation one arrives at 

= Co/3 {Bf + V{q) (f,^)) (5.3) 

which can be solved for the mean-field susceptibility x = ^it vanishing external 
field: 

^We^ss^^ = lim M = (5.4) 



The fluctuation-dissipation theorem (see e.g. ||56[, p. 116) gives an expression for 
the mean fluctuation of the order parameter: 

So far, no distinction between the correlated and uncorrelated contributions 
to the mean-fleld has been made; this we will now attempt, within the framework 
of Onsager's reaction-field theory. First, in the language of inflnitesimal 



variations of both the external and, thus, also the full mean field (starting from 
zero), the Weiss equation ( ^.l]) can be written as 

^^Weiss ^ ^^ext ^ J2 Vij5 ill,) 
j 

The Onsager prescription is to assume that the full mean field separates into 
two independent contributions, the correlated and the uncorrelated one, called 
reaction field and cavity field, respectively. The latter can be formally derived by 
subtracting the reaction field from the full mean field: 

SBr'' = + E V^^^ (/^.) - E ^^^^ (f^^),. ' (5-7) 
j j 

where the index "/ij" symbolizes the assumption that it is this part of the con- 
tribution to the field acting on /ij which is itself determined by the correlation 
between spin i and the surrounding spins j that are summed over. In a heuristic 
way, the value of 6 (fJ'j) ^. can, in principle, be determined from a kind of master 
equation written in terms of conditional averages: 

where {fj,j).^ ^ are the conditional average values of /ij for spin j when = ±1, 
and Piini) is the probability for spin i to have the magnetic moment /ij. Thus, 
since we are still dealing with the unmagnetized phase, we have 

p,(/i,) = ^(l±(/i.)). (5.9) 

The conditional averages are fixed values and so in ( p. 81) only the probabilities 
Pj(±l) can have a finite variation: 

= {6{fi,))M, (5.10) 

where in the last step we have again made use of the fact that we are dealing 
with the unmagnetized phase. Substituting this expression in the formula for the 
cavity field (|5.?1 ) gives 



(5.6) 



where we have introduced the parameter A as a measure of the interaction energy 
per spin: 



Y^V^.M- (5.12) 



By comparison with the earher Weiss result, we see that the reaction field has 
been isolated from the other contributions to the mean field, i.e. ^j^^eacuon _ 
+\5 {fii) , which must be subtracted as in when calculating the response of 



spin i to the Onsager cavity field. By way of Fourier transformation one gets 

5 (/i,-) = coP {5Bf + - \5 {fi^)) (5.13) 

and from this the susceptibility 

= ^_M = "A (5 14) 

In the light of the Weiss calculation this result could easily have been guessed; at 
this point, however, the calculation served the simple purpose of introducing the 
basic concepts of Onsager's reaction-field theory. 

5.1.2. The Heisenberg model 

We shall now apply the Onsager method to the Heisenberg model. While this 
generalization is quite straightforward, we want to point out one of the main 
shortcomings of the method, namely the fact that Onsager reaction-field theory 
is incapable of incorporating the effect of dimensionality on long-range order in 
a consistent way, contrary to what has been conjectured e.g. in or |]5^. We 



shall follow |5^ in treating the spin-1/2 anisotropic Heisenberg model given by 
^ = - E (^^^- (^r^J + ^'^D + ^^-^^^1) ' (5.15) 



n.n.ij 



where J^^ = eJy > Jij is assumed. This kind of anisotropy in spin space favour- 
ing uniaxial ordering is widely used (for a Green's function approach to be dis- 
cussed later within the Tyablikov approximation, see e.g. |31[|) and should, even in 



2D, certainly break the symmetry one encounters for the corresponding isotropic 
model. Once we have established the Onsager reaction-field theory for the Heisen- 
berg model ( |5.15D , we will be able to check the validity of the Onsager method by 
comparing its predictions with this expectation. 



Again, we start with a simple Weiss mean-field theory of an infinitesimal 
external field giving rise to an total mean field 

BY/''' = BtT + Y.J^3{S3)^ (5.16) 

j 

where the expectation value (Sj) is to be calculated self-consistently from 

^^i"^ ^ ^f3miB^<="'' (5-17) 

{S} 

and possible proportionality constants have been dropped for clarity (m: quantum 
number characterizing the z-component of the spin). Taking the condition of in- 
finitesimal external field seriously, we can expand the numerator and denominator 
of ( |5.17D to first order. 



^ m, (1 + (3miB^^''') 

{S} 

J2 (l + /3mi5]J--) 

{S} 

j:pf^yB^-, (5.18) 

where in the second step the denominator has been expanded to first order, i.e. 

-1 / \ -1 



5^ (1 + Pm^B^--) ^ 5^ 1 
{s} J \{S} 



1 



25 + 1' 



and we have also made use of ^ = 0. The sum over spin configurations, 

{s} 

symbohzed by '^^^ written explicitly as the sum over integer steps rrii = 

-S,-S + 1,...,+S : 



\m=—S / 

= {{m,)Y^' (3B^^^^^ (5.19) 



making use of the familiar partial sum formula ^^^i k"^ = and the 

isotropic distribution of spin orientations in the free case. Defining 

/=/3(Kf>=/3^^^^, (5.20) 



and using (|5.19|) and ( |5.16|) , we get for the Fourier transform of the spin z- 
component the result 



{SM) = 

leading to the Weiss mean-field susceptibility for the anisotropic Heisenberg model 

from which we can deduce an expression for the transition temperature by setting 
the denominator of (|5.22|) to zero so that the longitudinal susceptibility diverges, 



indicating a phase transition at the ordering wave vector q = Q: 

Tr = ^^J-\Q)- (5.23) 

Now, going over to Onsager react ion- field theory we write, by analogy 
with the Ising model, 

i i 

the second sum being again the reaction field which must be subtracted from the 
total mean field. The parameters \ij will in general be temperature-dependent 
and are determined by the correlation of the surrounding spins j with the "frozen" 
spin i. One easily derives the Onsager susceptibility 

xhnsi^ = i.^^i^j^i^^.xy (5.25) 



where A = J2j ^ij^ij- Assuming, as before, a phase transition at an ordering wave 
vector q = Q one gets for the critical temperature 

T^-' = ^^^^{J^{Q)-\) (5.26) 

which, dividing by ( p.23[ ), gives 

Singh 1^ has noted that a simple procedure to determine the parameter A 
appearing in ( |5.27| ) is provided by the fluctuation-dissipation theorem which, 
for zero external field, is 

/3 {SlSff = J2 X"{q)e''^'^^'~^'^. (5.28) 

In the isotropic case, one would have {{niif^ = ^^^^^^ and all values of rrii 
would have equal probability. In the anisotropic case larger values of \{Si)\ are 
favoured, while the corresponding pairs (S'f ) = ±m remain equally probable, so 
that (^{mif^ = a^^^^. Theorem ( |5.28| ) gives 



^ (5.29) 

With (ODD and (OD one finds x°(^f "') ■ J'iQ) = l£ which together with 
\ = JHQ)- JHQ) r^M F gives the final result for the Onsager susceptibility: 



_ TMF J.(g-) 



.0 



ax 



- TMF J.(g^ \ • (5-30) 



The constant x° cancels and, since we are in the spin-1/2 case, the constant a, 
which in the anisotropic case would normally "shift probability" to pairs (Sf) = 
±m with higher \m\, is unity, as=i/2 = 1, because ±1/2 are the only orientations 
possible. Thus, we get the final result 

'Y'Ons 1 
c 



TMF ^ 



Jf{9) ■ 



From this, we see that while Onsager reaction-field theory does modify the or- 
dinary mean- field result, it is independent of the anisotropy assumed above by 
setting = eJij, since the anisotropy parameter e will drop out of the frac- 
tion J^iq) / J^{Q) in the denominator. Effects of dimensionality that stem from 
the summation are, however, retained, e.g. its logarithmic divergence in two 
dimensions. 

In the next section, we shall discuss the same Hamiltonian using the Tyablikov 
method and find that it correctly treats this and other limiting cases. The at- 
tempt to enforce correct Mermin- Wagner type behaviour within mean-field theory 
obviously fails, at least as far as the Onsager improvement of mean-field theory is 
concerned. We note that this observation refutes contrary assertions in the 



literature (such as []57|). 



5.2. Mermin- Wagner Theorem and Tyablikov Method 

Tyablikov proposed a method of approximation for the Heisenberg model by 



suggesting a decoupling in the higher-order Green's functions derived within the 
equation of motion method. Let us remind ourselves of the Heisenberg Hamilto- 



H = -Y, MStSj + StS^^) - gj^B, J2 St, (5-31) 

ij i 

where we have added a symmetry-breaking term Hf, = Qj^Bo Sf := b S'f to 
the unperturbed Hamiltonian Hq. The coupling constants Jij fulfill the symmetry 
conditions {Jij = Jji] Ja = 0) and are taken to give the ferromagnetic case: 
Jij > 0. 

In order to evaluate the problem, we employ the equation of motion method 
by use of Green's functions. Since we are interested in the question whether or not 



the system described by the many-body model Hamiltonian exhibits spontaneous 
ferromagnetic order, we choose as the appropriate Green's function 

Gif{t-t') = {{st{ty,s-{t'))r' 

= -ze{t-t'){[Snt),S.{t')]_) (5.32) 
which, when plugged into the first-order equation of motion gives 

= 2h% {St) + hhGf{E) - 2hY,J^m {{{s^st- s-)y;' 

m 

-{{s'r.sr;Si)y;'). (5.33) 



The higher-order Green's functions must be decoupled in order to make the 
problem analytically and numerically tractable. In the Tyablikov decoupling 
scheme, this is done by separating from the spin-flip operators in the 
higher-order Green's function ad hoc via the substitution 

{{S^^Sr, S-)) (5^) {{Sr, S-)) . (5.34) 

One then arrives at the equation 

(e - gjtXBB, - 2h {S^) J2 Jrm) Glf{E) = 2h%, {S^ - 2hY,{St) GZ]{E)J,^ 

\ m / m 

(5.35) 

which can be solved for the Green's function explicitly in momentum space with 
the standard Fourier transforms {q is chosen from the first Brillouin zone) 



1-3 



Jij 



(5.36) 
(5.37) 



thus arriving at the solution 



where the energy 

E{q} = 2h {S') (Jo - J{q}) + gjfisBo (5.39) 

has been introduced {Jq := J2i Jij)- Note that in the transition from eqn. ( ^.35|) 
to eqn. (|5.38|) we have made use of translational symmetry 

{St) = (5^,) = (S^) , (5.40) 

thus assuming the ferromagnetic case where all spins point into the same direction. 
The Green's function, by virtue of the relation S(^{E) = —^JG^^^^{E) gives the 
spectral density 

S,iE) = 2h^ {S') S{E-2h {S') (Jo - J{q)) - gjfisBo) (5.41) 

which allows one to make use of the spectral theorem for calculating the expecta- 
tion value 

(^7S?-> = i E £ dE^. (5.42) 



q 



In the diagonal case {i = j), spin algebra yields (restricting our attention, for the 

moment, to the S = ^ case) the formula (S'~S'.^) = h^S — hl^S^^-^ ■ (Because 

of translational symmetry ( |5.40|) we shall drop the site index i.) Combining the 
two formulas, a compact though implicit equation for the expectation value of the 
2;— component of the spin is derived: 

Up to this point, we have only summarized the well-known Tyablikov procedure for 
calculating the magnetization. In the following section, we shall examine whether 
or not the Tyablikov decoupling conserves the Mermin- Wagner theorem, which 
states that there can be no finite-temperature (/? < 00) magnetic phase transition 
in one or two dimensions. 



5.2.1. Tyablikov decoupling and Mermin- Wagner theorem 



The standard proof of the Mermin- Wagner theorem for a given model Hamilto- 
nian starts with the Bogoliubov inequality, from which an upper bound for the 
magnetization is derived which is then shown to vanish in the limit i?o ^ at 
any finite temperature. In our case, discussing an approximative method, the Bo- 
goliubov inequality cannot be used for testing the Mermin- Wagner theorem since 
the method does not approximate the Hamiltonian but the higher-order Green's 
functions instead. However, since with eqn. ( p. 431 ) an expression for the magneti- 



zation (5*^) is given, it is only natural to examine this equation in order to check 
whether the Mermin- Wagner theorem remains indeed valid within the Tyablikov 
method. 



First we reformulate (|5.43|) writing it as 



Tj E _ .-m^n =irjY. '-E{q) (5.44) 



g g 

which, in the continuum notation appropriate for the thermodynamic limit, can 
be written as 

^^'^-77^ / d''qcoth(^^i2h{S^)iJo-J{q))+gjfiBBo)y (5.45) 

l.B.Z. 



where the full expression for the energy E{q) has been inserted and v^'^'' = lim 

is the specific (c?— dimensional) volume. A possible phase transition can reasonably 
be expected only for temperatures below a certain critical value T^, for which the 
limit of a vanishing external field must then be considered: 

5o = 0+;T-T, = 0" ^ (5") = 0+. (5.46) 

Using the series expansion of the hyperbolic cotangent in the assumed region of 
finite spontaneous magnetization, 

oo 

cothx = 1 + 2 exp (— 2mx) , (5-47) 

m=l 



one finally arrives at 

(S^ ^ 1 

1 + 2 ^ e-t^^^aji^BBoJ^ J rfrfg-exp {-2mph {S') (Jo - J{q))) 

m=l ^ l.B.Z. 

(5.4^ 



1 + 2 E 



-l3mgjfj,BBo . j(d) ^ + 20^''^ 



m=l 



(which at the same time defines the quantities I^^^ and (p^'^^). 

The choice of dimension D will turn out to be crucial for evaluating the be- 
haviour of (S^) in the limit Bq ^ 0. 

5.2.2. Proof of Mermin- Wagner theorem in two dimensions 

To prove the Mermin- Wagner theorem, it must be shown that in two dimensions 
there exists a bound to (S^) that vanishes in the limit Bq ^ (see eqn. ( [5 .461 )). 
This is certainly the case if the quantity 

°° (2) r 

0(2) = ye-^'"^'^^^^"-^ / d^qexp{-2mPh{S')iJo-Jm) (5.49) 

l.B.Z. 

diverges. 

In a first step, the support of the integral is restricted to a sphere inscribed 
into the first Brillouin zone with radius go- Secondly, the integrand itself can be 
bounded from below by making use of the fact that Jq — J{q) < Qq^, as long as 

the interaction decays rapidly enough for the quantity Q = ^ Yliij ~ l^yl 
to converge. Now spherical coordinates can be used, and since the integrand has 



been replaced by a function of \q\ only, one finds = const.) 



,(2) 



10 



m=l 
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(-1) 







,(2) 




m 



m 



(5.50) 



Assuming for the moment a finite value of (5*^) , the second sum converges in 
the limit -Bq — >■ 0, since subsequent terms vanish exponentially with m. The first 
sum, however, is the harmonic series and, thus, diverges. But as 0'-^'' diverges, 
eqn. ( ^.48|) states that in this case {S^) must be zero. Starting from a finite value 
of {S^) one arrives at a contradiction, so the assumption must be wrong. 

At first sight, it may appear obvious from (|5.5CI|) and (|5.48|) that (5*^) = 
is indeed consistent with both equations, since the two series appearing in the 
lower bound of 0*^^-' cancel exactly. However, in the derivation it has been tacitly 
assumed that (5*^) ^ 0. This can be seen by explicitly inserting (|5.5CI| ) into ( |5.48D 
to get the full relation for the magnetization: 



hs 



< 



1 



1 + 



47r/3?iQ(5^> 



^ exp {-(3mgjfiBBo) 1 - exp -2m(3hQ (5^) 

n=l L \ /J 



(5.51) 



where the assumption {S^) = would lead to an expression "0/0" in the denom- 
inator. Thus, this case can only be analyzed by means of I'Hospital's rule, where 



one considers the limit lim(5z)^o liniBo^o ("■■)• 



lim lim 



< 



1 



^1 {s->-o 



(5.52) 



where the last equality is due to I'Hospital's rule. Reformulating the initial equa- 
tion in this way, the denominator now is explicitly divergent, so indeed (S^) = 0. 

The situation changes when an external field i?o 7^ is applied: As one can 
easily see from ( |5.50| ), both series then converge; thus is also convergent and one 
arrives at a non-vanishing upper bound to the magnetization. The upper bound is 
still a function of (S^) , so it will be difficult to extract more detailed information 
about (S^) from this relation except for the mere possibility of a finite value of 
the magnetization. 



5.2.3. Three-dimensional case 

Testing the above procedure in the 3-dimensional case provides another valuable 
check. We shall, at this point, show that the lower bound of derived by analogy 
with the two-dimensional case does not diverge in three dimensions. This is, of 
course, only a negative statement about the particular lower bound. Typically, for 
tests of the Mermin- Wagner theorem in three dimensions, this is all one can expect, 
since a positive statement about a possible lower bound of the magnetization (i.e. 
a converging upper bound of 0) will rarely be possible. Later, we will present 
a different treatment of the three-dimensional case which will also allow us to 
discuss the dimensional crossover from two to three dimensions (or vice versa) in 
much greater detail. For the moment we note that, indeed, the finite upper bound 
to the magnetization survives the limit i?o — ^ in three dimensions: 



hs 



m=l ^ ^ 

(°° 7,(3) y / . N 
1 + 2 ^ e-^"^^^''^^"— / dqq^ exp (^-2m/3/i (5") Qq^) 
m=l ^ 



V 



due-"" 



(3) 



= |^l + 2^^e-'^'"^^'*«^°^A^j (5.53) 
with the quantity defined as 



-1 f d^ \ d 



— / due-"" 



2 



/ \ ^/2l3mhQ{S-)qo 

d 1 



where we have made use of the chain rule for the new variable 

e = ^J 2PmhQ {S^)qo. 



J due-""', (5.54) 



X 

The first integral is of the form ^ J f{y)dy = f{x), so that 

^ 







due " 



(5.55) 



Let us discuss this equation in relation to the inequality ( ^.53| ) . Assuming a finite 
non-vanishing value of the magnetization, we choose a fixed non-zero number e 
such that (5*^) > e > 0. We can then approximate the second integral in ( |5.55| ) 



from below by setting the integrand equal to the value of e " at its upper bound, 
thus strengthening the inequality in (|5.53| ): 



Am > ^ (e - 1) ^ ^ (5.56) 

The series in the denominator of the R.H.S. in ( ^.53| ) then always converges, even 
in the limit Bq 0, so that a magnetic phase transition is not ruled out. This is, 
of course, no proof that a magnetic phase transition must occur; as one can easily 
see from ( 5.55| ) by taking the limit (S^) 0, the individual terms of the series 



vanish: —— 0, which is not in contradiction with ( |5.53| ). 

(d ) ^ 

5.2.4. Absence of spontaneous sublattice magnetization in ABAB-type 
antiferromagnets 

While the Tyablikov method has been very popular for treating ferromagnets, 
where it has led to convincing results over a broad temperature range, it can, in 
fact, be generalized to cover antiferromagnets as well [Q, |^0[. Recently, there has 



been renewed interest in the antiferromagnetic Heisenberg model with relation to 
high-temperature superconductors such as YBa2Cu306-x which, for low doping, 
displays antiferromagnetic behaviour. 

Quite generally, in the Heisenberg antiferromagnet, the respective sublattices 
can be treated using the Tyablikov formalism essentially in the way discussed 
above for ferromagnets. As to antiferromagnets, we want to sketch briefiy how 



to proceed. Following Hewson and ter Haar pO| we consider, for simplicity, the 



case of the lattice being divided into two sublattices A and B so that the near- 
est neighbours of every atom of one sublattice all belong to the other sublattice 
(ABAB-type ordering). (More complicated types of antiferromagnetic order can 
be discussed in a similar way; for layered antiferromagnets see e.g. ref. [^.) If 



we restrict our attention to the case of nearest-neighbour interaction, thus keep- 
ing only the essential features of the interaction, the Hamiltonian (|5.31|) can be 
written as 

« = + S}S;) - (j2 SI - S') (5.57) 

f g \ g f / 

with / and g referring to the sublattice with predominantly spin-down or spin-up 
orientation, respectively. (The sign change for |-spins coupling to the external 
field is a matter of convention.) 

Since we are now dealing with antiferromagnetic order, we can no longer as- 
sume translational symmetry as in eqn. ( |5.40| ). Instead, for purely geometric rea- 



sons, we must distinguish between "same-lattice" Green's functions (^(^S^,Sj,')') 
or ((5'+; 5^)) (where /, /' and g,g' denote sites of the same sublattice A or 
B) on the one hand and "mixed" Green's functions (^(^S'^; Sj)) or ((S'j;S'^)) 
on the other hand. The equation of motion ( p.33| ) must be modified in order to 



accommodate the distinction between the two kinds of Green's functions. Thus, 
one arrives at a system of two equations of motion (one for the "same-lattice" 
Green's function, the other for the "mixed" Green's function) which can be writ- 
ten in matrix form. In the "mixed" case there will be no inhomogenity as spin 
operators at different lattice sites commute (which will always be the case when 
dealing with different sublattices). The Fourier transforms of "same-lattice" and 
"mixed" Green's functions (denoted by G^j^ and G^j:) are taken over the respective 
reciprocal sublattices in fc-space: 

Gr, = j:((S:-^S,))e-^H^^-^.') (5.58) 
g~g' 

= (5-59) 

g-f 

where the summation is over sites of the same sublattice or over different sub- 
lattices, respectively. The higher-order Green's functions are decoupled, as in the 
ferromagnetic case, through the Tyablikov procedure ( |5.34D . Inverting the 2x2 
Green's function matrix, one can solve for the same-lattice Green's function G^j^, 



which, via the spectral theorem and the identity (5*^ S^) = h^S — h (^S^-^^ for the 
S = ^ case, gives an expression for the sublattice magnetization: 



hS 



(5.60) 



The formula resembles closely that for the ferromagnet, except for the typical 
square-root dependence on the exchange integrals which is related to the spin-wave 
dispersion in simple antiferromagnetic spin- wave theory 0. (For more details on 
the derivation see also refs. and [^].) Earlier, for reasons of simplicity, we 



had restricted the interaction to nearest neighbours. Thus, 

n.n. 

i=fix 

Jio) = 5^ V^"^"^'-"^^), (5.62) 

n.n. 
i=fix 

where ^^n.n." stands for "nearest neighbours only" and i =fix reflects the fact that 
only the dummy variable j is summed over. 

From ( [5.60[ ) it is possible to derive an explicit formula for the RPA Neel tem- 



perature T§^^, by making use of the series expansion coth(x) ^ + ... for small 
values of x, i.e. close to the transition point where {S^) (T = T/v) = : 



Jn _ 1 

2 



= g^ 2v^ ^ 1 • (5.63) 



k 



In our case, where Jij = J 7^ for nearest neighbours i, j only and where the lattice 
is of ABAB structure, it is straightforward to show that in the thermodynamic 
limit the sum in the denominator diverges in two dimensions while remaining 
finite in three dimensions. Due to a l/k"^ divergence of the integrand close to 
k = 0, which is not cancelled by the 2D volume element kdk, the sum diverges 
logarithmically in two dimensions. In three dimensions, the sum remains finite. 
Thus, one recovers the Mermin- Wagner theorem which states that there can be 
no finite-temperature phase transition in the two-dimensional isotropic Heisenberg 
model. 



5.2.5. Transition from two to three dimensions via 3D anisotropic in- 
terlayer coupling 



Previously, we saw (eqn. (^.48|) , (|5.6CI|) ) that the dimension determines the (sub- 



lattice) magnetization via summation over £)-dimensional fc-space. In those cases 
where one starts from an a priori restriction in geometry, the choice of dimension 
is fixed from the very beginning. These cases were discussed at length for the 
ferromagnet and, in the previous section, for the antiferromagnet as well. There 
is, however, another method for changing effective dimensionality. Instead of re- 
stricting the geometry of the system a priori, one can change the interaction in 
such a way as to mimic the change in dimension. 

As an example, we shall return to our discussion of the antiferromagnet. The 
interaction had been restricted to nearest neighbours. If we choose the interaction 
parameters Jij so that (for a fixed lattice site i) Jij = J|| if sites i and j lie within 
a plane and Jij = J± if j is located in the direction orthogonal to this plane, then, 
by varying the anisotropy parameter 



J± 
J\\ 



(5.64) 



we can change the coupling continuously from quasi two-dimensional {e = 0) to 
(isotropic) three-dimensional {e = 1). 

Of course we have to distinguish between the coordination numbers within the 
plane and those along the orthogonal direction, z\\ and zj_. Modifying ( 5.63| ), the 



formula for the RPA Neel temperature, to take into account the anisotropy gives 
(see also ref. |^ for comparison) 

tRpa _ J\\ zii+ez^_ 

~ 2k~ ^ V ^ [0.00) 

where the substitution Jq = Yl Jij ~ + ^^±J\\ has been carried out and 

n,m 
i=fix 

where we have used the common abbreviations 

711 (g) = f E^^''^" (5-66) 

II n.n.\\ 

^-L I 
n.n.± 



Note that in the strictly two-dimensional case (where z± = Oand q = {qxily)) 
the sum in the denominator ( ^.651 ) diverges as expected. In the following, however, 
we shall focus on the limit with _D = 3, i.e. q = {qx,qy,qz) = (^|,'?±),but with 
quasi two-dimensional interaction: £ — >■ 0. 

Apart from prefactors, in the thermodynamic limit the sum can be replaced 
by an integral: 

d^q ; ^. (5.68) 

^ z||7||(g)+£z±7±(g) ^ 

2||+£2j_ 

Evidently, the main contribution to the integral will come from near q = 0, since 
there 7||(g) and 7±(g) can be approximated by 1 — a\\/±q'^. Note that by a proper 
choice of a||and a_Lone can isolate a finite sphere around the origin throughout 
which the approximation (as far as the total value of the integral is concerned) 
holds exactly. With a suitable choice of z\\, it is also possible to write 

II s 

n.n.H 

instead of the original definition, while preserving the qualitative dependence on 
\q\ (at least in the thermodynamic limit.) 

With the above modifications, we find for the integral /over a suitable sphere 
around the origin 



s 



+ 4(^)'(l-a||,^)(l-«,,i)^ 



where only quadratic contributions have been retained in the last line. Quantities 
with e as an index are smooth functions of e remaining finite as e — ^ 0. We have 
also made extensive use of our previous notation q= {q\\,q±) where q\\ is a vector 
of the basal plane and q± the coordinate orthogonal to the plane. 



The evaluation of the last integral is carried out in cylinder coordinates, first 
integrating in ^-direction from a lower cutoff — Aj^ to the upper cutoff +A±. Af- 
terwards, we integrate in the basal plane using polar coordinates: 



I' 



(fq / dq 



1 



(^A + B,{1 - f3e^)ql^ + e^C^ql 



"II 

2n J dq\\ 



2?|| 



arctan 



eVCWAj 



Ae + B,{l-f3e^)ql 



IA, + B,{l~f3e^)ql 



A-k/ [e^/Ce^K^ f ^ arctan M 



{5.71] 



Uo 



where in the last integral the substitution 



u :- 



A, + 5, (1 - I3e^) q^ 



has been used, with mq = oo and ui 



M,+B,(l-/3e2)A2 



(5.72) 



as the integration limits. 



The integral vanishes at infinity and, thus, we find 
47r/ [eVCTi^A.) 



I' 



A, + B,{l-l3e^)K^ 



arctan 



ev^Aj 



1 / A, + i?,(l-/?£2)A 
+ o In U + 97T—^ 



lA, + B,{l-l3e^)A^^ 

(5.73) 



From this, one can read off the following result: While the arctan approaches zero 
for £ — s> and thus may cancel the divergency due to the prefactors, there 
will always remain a divergency in the logarithmic contribution. Since the integral 
/' determines the Neel temperature via eqn. (|5.65| ) 



rpRPA 

-'at 



J\\{z\\+ezA_] 
ikB 



(5.74) 



we see that in the quasi-2D hmit of vanishing interlayer couphng [e —>■ 0) the Neel 
temperature (as calculated from Tyablikov RPA) vanishes as 

in complete accordance with the Mermin- Wagner theorem. 

5.2.6. Extending the results to higher spins (S > 1/2) 

The results so far were derived for the 5 = 1/2 case only. They can, however, be 
extended to arbitrary spin by a rather simple procedure. The definitions of 0^'^^ 
and E{q) remain unchanged in relation to the previous discussion for ferromagnets 
|]6^ as well as for antiferromagnets Only eqn. ( ^.48[) must be modified to 



run 

{s^ ^ [s - HS)] [I + HS)r^' + {1 + S + <P{S)) {HS)f'^' 
h (1 + 0(5))'^+' - (0(5))'^+' 

The easiest way to reproduce the above 5 = 1/2 results for higher spins is to show 
that the convergence, or divergence, of 0^'^^ will lead to a finite or vanishing {S^) 
value, respectively. This suffices because 0^'^^ is the only quantity that depends 
explicitly on the dimension. 

For large values of 0^*^^ one can expand ( [5.76| ) in powers of 1/0(^*^0 get 

{S^) ^ S{S+l) 



thus establishing to order (9(1/0^) that the results obtained above carry over to 
this section. 



Alternatively, writing ( p.76| ) in powers of 0, one finds 

js^) ^ [s - 0] [1 + (pr^' + (1 + g + 0) 0^^+^ 

h ~ (1 + 0)2^+^ - 025+1 

2S-\-l 25+1 

_ n=0 n=0 

~ 2S+1 

n=0 

25 



n=l 



25 
n=0 

25—1 
n=l 



25 



(5.78) 



n=0 



where in the last step we have used the recursion relation for binomial coefficients. 
For S > 1 the coefficients in the numerator and the denominator are finite, so 
from the general form 



25—1 
E «n0" - S 

(5.79) 



h 25 

we can see that the polynomial in the denominator is of higher order than the nu- 
merator since (32S = (^25^) 7^ ^- ^ diverges assuming finite {S^) , eqn. 
( ^.79| ) immediately contradicts the assumption. Thus, in two dimensions, (S^) 
must vanish. On the other hand, if remains finite (as expected in 3D), a finite 
value of (5*^) is consistent with ( |5.79| ), thus allowing for a non-zero spontaneous 
magnetization. 



To summarize our discussion of the Tyablikov method, it has been shown that 
the Tyablikov method for the Heisenberg model conserves the Mermin- Wagner 



theorem: In two dimensions, a finite value of (S^) can only be achieved by apply- 
ing an external field; a magnetic phase transition at finite temperature is ruled 
out. Spatial dimensionality enters the calculation explicitly through the fc-space 
integral appearing in the definition of (p. The Tyablikov method, in this respect, 
differs markedly from mean-field theory, where the coordination number of a given 
lattice is the only variable accounting for the geometry of the system, or from sim- 
ple improvements thereof, such as Onsager reaction-field theory discussed in the 
previous section. 



6. Summary 



In this paper, we have investigated several aspects of dimensionahty in many- 
body theories of long-range order. The Bogoliubov inequality, together with the 
concept of quasi-averages, which we have discussed in Chapter II, provides the 
mathematical tools for examining rigorously whether or not a phase transition 
occurs. For film systems, i.e. systems of d identical two-dimensional layers stacked 
on top of each other, we have shown that a magnetic phase transition at finite 
temperatures can be rigorously ruled out, unless d — > oo. The proof holds for a 
wide class of many-body Hamiltonians, such as the Heisenberg, Hubbard, Kondo- 
lattice, and Periodic Anderson Model. The general nature of the proof suggests 
that analogous results would be obtained for other Hamiltonians and other kinds 
of phase transition as well. For s-wave pairing in Hubbard films this conjecture 
has been confirmed. 

The Bogoliubov inequality has been adapted to fit in with a more general 
framework of correlation inequalities as presented in Chapter IV. From the alge- 
braic properties of some of the main statistical quantities characterizing a phase 
transition (most notably the static susceptibihty) , rigorous relations have been 
derived which may be used to derive relations between different types of order 
parameters or phase transitions. The general applicability of the formalism pre- 
sented, particularly the fact that the Bogoliubov inequality and the Coldstone 
theorem fit in well with it, suggest that further research in this direction might 
be able to clarify some conceptual issues in the theory of (magnetic) phase tran- 
sitions. 

Finally (Chapter V), we have contrasted some approximative methods of 
many-body physics with respect to whether or not they conserve the Mermin- 
Wagner theorem. As it turns out, simple extensions of mean-field theory (which 
is known to violate the Mermin- Wagner theorem) satisfy the Mermin- Wagner the- 
orem at best on a superficial level. By including fiuctuations as in the Tyablikov 
method, it is possible to construct an approximative theory that reproduces the 
Mermin- Wagner theorem in all the relevant limiting cases. 
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As of today, the Mermin- Wagner theorem is perhaps the strongest rigorous 
statement available, as it covers a variety of different models, is independent of 
the microscopic parameters and is applicable in the thermodynamic limit. It is, 
however, a negative statement about phase transitions. Therefore, it is unclear 
whether or not future investigations of this problem will eventually be able to 
satisfactorily resolve the question of precisely how phase transitions do occur in 
the real world. With these considerations in mind, it is perhaps appropriate to 
close our discussion of the Mermin- Wagner theorem and related exact results on 
a philosopher's remark, namely that 

"It is the mark of an educated man to look for exactness in each class 
of things just so far as the nature of the subject admits. "[| 



Aristotle, Nicomachean Ethics. 



Bibliography 



[1] Redhead, M.: Brit. J. Phil. Sci. 31 (1980) 154 

[2] Nolting, W.: Quantentheorie des Magnetismus, Vol. 1, Teubner: Stuttgart 
1986. 

[3] Nolting, W.: Quantentheorie des Magnetismus, Vol. 2, Teubner: Stuttgart 
1986. 

[4] Domb, C: The Critical Point, Taylor & Francis: London 1996. 

[5] Heisenberg, W.: Z. Phys. 38 (1926) 441 

[6] Hubbard, J.: Proc. Roy. Soc. A267 (1963) 238 

[7] Gutzwiller, M.C.: Phys. Rev. Lett 10 (1963) 159 

[8] Kanamori, J.: Progr. Theor. Phys. (Kyoto) 30 (1963) 275 

[9] Nolting, W.: Grundkurs: Theoretische Physik, Vol. 7 (Viel-Teilchen- 
Theorie), Vieweg: Braunschweig/ Wiesbaden ^1997. 

[10] Bogoliubov, N.N.: Physica (Suppl.) 62 (1960) S 1 

[11] Bogoliubov, N.N.: Phys. Abhandl. Sowjetunion 6 (1962), 1, 113, 229 

[12] Akhiezer, A.I. and S.V. Peletminskii: Methods of Statistical Physics, Pcrga- 
mon Press: Oxford 1981. 

[13] Hohenberg, P.C.: Phys. Rev. 158 (1967) 383 

[14] Mermin, N.D. and H. Wagner: Phys. Rev. Lett. 17 (1966) 1133 

[15] Mermin, N.D.: Phys. Rev. 176 (1968) 250 

[16] Fernandez, J.F.: Phys. Rev. A2 (1970) 2555 

89 



[17] Kishore, R. and D. Sherrington: Phys. Lett. A 42 (1972) 205 

[18] Chester, G.V.; Fisher, M.E.; and N.D. Mermin: Phys. Rev. 185 (1969) 760 

[19] Wegner, F.: Phys. Lett. 24A (1967) 131 

[20] Walker, M.B. and T.W. Ruijgrok: Phys. Rev. 171 (1968) 170 

[21] Ghosh, D.K.: Phys. Rev. Lett. 27 (1971) 1584 

[22] van den Bergh, M. and G. Vertogen: Phys. Lett. A 50 (1974) 85 

[23] Robaszkiewicz, S. and R. Micnas: phys. stat. sol. (b) 73 (1976) K35 

[24] Baryakhtar, V.G. and D.A. Yablonskii: phys. stat. sol. (b) 72 (1975) K95 

[25] Thorpe, M.F.: J. Appl. Phys. 42 (1971) 1410 

[26] Krzemihski, S.: phys. stat. sol. (b) 74 (1976) K119 

[27] Matayoshi, G. and S. Matayoshi: Bull. GoU. Sci. (University of the Ryukyus) 
64 (1997) 1 

[28] Rastelh, E. and A. Tassi: Phys. Rev. B 40 (1989) 5282 

[29] Li, Y. and K. Baberschke: Phys. Rev. Lett. 68 (1992) 1208 

[30] Farle, M.; Baberschke, K.; Stetter, U.; Aspelmeier, A.; and F. Gcrhardter: 
Phys. Rev. B 47 (1993) 11 572 

[31] Haubenrcisscr, W.; Brodkorb, W.; Gorciovei, A.; and G. Gostache: phys. 
stat. sol. (b) 53 (1972) 9 

[32] Wu, J.H. and W. Nolting: phys. stat. sol. (b) 219 (2000) 181 

[33] Gelfert, A. and W. Nolting: phys. stat. sol. (b) (2000) 805-818 

[34] Gelfert, A. and W. Nolting: J. Phys. Gond. Mat. (in press) 

[35] Proetto, G. and A. Lopez: J. Physique Lettres 44 (1983) L635-L640 

[36] Gassi, D.: Physica A 191 (1992) 549 

[37] Uhrig, G.S.: Phys. Rev. B 45 (1992) 4738 



[38] Noce, C. and M. Cuoco: Phys. Rev. B 59 (1999) 7409 
[39] Yang, C.N.: Phys. Rev. Lett. 63 (1989) 2144 

[40] Su, G.; Schadschneider, A.; and J. Zittartz: Phys. Lett. A 230 (1997) 99 

[41] Chandra, P. and P. Coleman: Int. J. Mod. Phys. B 3 (1989) 1729 

[42] Shraiman, B.L and E.D. Siggia: Phys. Rev. Lett. 62 (1989) 1564 

[43] Hattori, K.; Hattori, T.; and H. Watanabe: Progr. Theor. Phys. (Kyoto) 92 
(1987) 108 

[44] Auerbach, A.: Interacting Electrons and Quantum Magnetism, Springer- 
Verlag: Berhn 1994. 

[45] Kubo, K. and T. Kishi: Phys. Rev. B 41 (1990) 4866 

[46] Pitaevskii, L. and S. Stringari: J. Low Temp. Phys. 85 (1991) 377 

[47] Lange, R.V.: Phys. Rev. 146 (1966) 301 

[48] Lange, R.V.: Phys. Rev. 156 (1967) 630 

[49] Brenig, W.: Z. Phys 169 (1962) 219 

[50] Wagner, H.: Z. Phys. 195 (1966) 273 

[51] Hohenberg, P.C. and W.F. Brinkman: Phys. Rev. B 10 (1974) 128 

[52] Stringari, S.: Phys. Rev. B 49 (1994) 6710 

[53] Onsager, L.: J. Am. Chem. Soc. 58 (1936) 1486 

[54] Brout, R. and H. Thomas: Physics 3 (1967) 317 

[55] Singh, A.: arXiv.org (2000): |cond-mat/OOlIMl| 

[56] White, R.M.: Quantum Theory of Magnetism, SprmgeT-Verlag: Berhn ^1983. 
[57] Tusch, M.A.; Szczech, Y.H.; and D.E. Logan: Phys. Rev. B 53 (1995) 5505 
[58] Bogohubov, N.N. and S.V. Tyabhkov: Soviet Phys. - Doklady 4 (1959) 604 



[59] Pu, F.-C: Soviet Phys. - Doklady 5 (1960) 128 

[60] Hewson, A.C. and D. ter Haar: Physics 30 (1964) 890 

[61] Lines, M.E.: Phys. Rev. 131 (1963), 540 

[62] Lines, M.E.: Phys. Rev. 135 (1964) A1336 

[63] Lee, K.H. and S.H. Liu: Phys. Rev. 159 (1967) 390 

[64] Praveczki, E.: Phys. Rev. 130 (1963) 147 

[65] Singh, R.P and M. Singh: phys. stat. sol. (b) 169 (1992), 571 

[66] Gradshteyn, I.S. and I.M. Ryzhik: Table of Integrals. Series, and Products, 
ed. A. Jeffrey, Academic Press: San Diego/London '''1994. 

[67] Hardy, G.H.; Littlcwood, J.E.; and G. Polya: Inequalities, Cambridge Uni- 
versity Press: Cambridge 1934. 

[68] Luttinger, J.M.: Progr. Theor. Phys. SuppL (Kyoto) 37/38 (1966) 35 



A. Positive semi-definiteness of the 

BOGOLIUBOV SCALAR PRODUCT 



With ( p.l3|) the Bogoliubov scalar product ?B has been defined as 

^{A-B)= J2 {n\A+\m) {rn\B\n) '"^ (A.l) 

En^Em 

with 

_ exp(-/3E„) 
'^"-tr(exp(-/3i/))' ^^-^^ 

where H is the Hamiltonian and A, B are arbitrary operators. 

In order to verify that 03 indeed has the properties one expects from the 
mathematical definition of a scalar product, we check the axioms: 

• 03(^4; i?) is a complex number with Q3(A;i?) = (?B(-B;A))*. This follows 
from ( |A.1| ) by noting that {wm — Wn)/{En — Em) is a real number and 

{^(n \B^\ m) (m \A\ n))* = (n \A'^\ m) {m \B\ n) . (A.3) 

• ^{A]B) is linear with respect to its arguments: ^{A^aiBi + 02-82) = 
ai^{A] Bi) + a2'*&{A\B2). This follows from the properties of the matrix 
element {m\B\n) . By use of the first axiom, one can derive an analogous 
relation for the first argument of *B. 

• 03 induces a positive semi-definite norm in the space of operators, i.e. 

Q3(A;A)>0, (A.4) 
because of {Wm — Wn)/{En — Em) > 0. 

• ^{A; A) = does not, in general, imply A = (for which (|A.4|) indeed 
vanishes); e.g. for the Hamiltonian, ^{H;H) = , while = is no 
physically meaningful Hamiltonian. 
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B. Upper bound for 



c 



+ 

nf3a 



In order to derive an upper bound for the expectation value |(c^/3o-Cfc^o-) | one uses 
the procedure that has been discussed in the main text (see eqns. ( p.39|) ff.) for 



the correlation function (c^^^c„^_o.5'„^) which we encountered in the s-f case. 
For (n, /3) = (^,7) one obviously has 

(c+^.c„^.> < 1. (B.l) 
Starting, in the general case, with the decomposition 

~ i^iaa ~ ^Iva) {'^iaa ~ Qua) 

-i {4aa - ictua) i^iaa + iClua)) (B.2) 

one can apply the spectral theorem, 

4 00 

—00 

where : (j = 1,2,3,4) — * (0(j) = +1; — 1, +^5 — defines a phase factor for 



every term as in (|3.39|) , and S + (E) = ^ ( [Aj^, Aja] ) (E) is the spectral 
density in its energy representation. Hence 

4 °o 



4 — 
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where each of the 0-th spectral moments in 
tation value of the anticommutator: 



can be expressed by the expec- 

—oo 

In ( [B.4| ), use has already been made of the triangle inequality 



(B.5) 



< 



and of the fact that for pairs of adjunct operators, the spectral density is positive 
definite. 

Since we have discussed the case {n, P) = (^,7) in ( |B.1D separately, we can 
now assume (n, /?) 7^ (A;, 7). The 0-th spectral moments can be easily calculated; 
they each give M^^\+ = 2. Thus for all n, k, (3, 7 we have as an upper bound 



\^n0(7^k~i(j) < 2. 



(B.6) 



C. From correlation inequalities to 

THE BOGOLIUBOV INEQUALITY 



For completeness, we derive in this section the remaining integrals in eqn. ( |4.72| ), 
thus completing the derivation of the Bogoliubov inequality within the correlation 
inequalities approach: 



For the integral on the LHS of ( [4.72| ) it must be shown that 



h I ~^^\H,C]_;[H,C]+SE) 



(C.l) 



Since the double commutator is related to the spectral density Scc+ by 



dEESccAE), 



(C.2) 



due to the general commutator representation of the spectral moments, it 
suffices to show directly that 



H,C]_;[H,C]: 



iE) = E^Scc+iE). 



(C.3) 



This is indeed the case: 



[H,C]_;[H,C]\ 



i J2 (e-^'''" - e-^^") (m \HC - CH\) ■ 



{n\{HCy - iCHy\ m) 5{E - (E„ - EJ) 

1 J2^e-^E^ _ e-/3^") - Ejf {m \C\n) ■ 



mn 



\C+\m)5{E-{En-E^)) 
i (e-^^" - e"^^™) E^ (m \C\ n) {n\C+\ m) 



■5{E - (E„ - E^)) 
E^Scc+{E), 



(C.4) 
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where we have twice made use of the spectral representation (|4.38| ) of the 
spectral density. 



The integral on the LHS of ( [4 .721) can be treated similarly 

/dE ^ 
-^Slc,H]_;A+{E) 

^ f^J2 (e-^'''" - e-^^") (m \CH - HC\ n) ■ 

mn 

{n\A+\m)6{E-{Er,-Ej)\^ 

4e( 



) (m \C\ n) {n\A+\m)5{E- (E„ - E^)) 



dEScA+{E) 



(C.5) 



Thus the transition from the correlation inequality ( [4.72| ) to the Bogoliubov 
inequality (^4.73|) can indeed be made. 



D. Integral inequalities 



We give here, for reference, a list of important mathematical inequalities that 
hold for real-valued functions /(x), g{x) under the additional conditions given 
with each inequality.^ Not all of these inequalities are used in this thesis, but it 
appears worthwhile to contrast those we do use with related ones. 

Cauchy-Schwarz inequality 

Let f{x) and g{x) be any two integrable functions on [a, b] . Then, 



dxf{x)g{x)j < dxfix)j dxg\x)j , (D.l) 
and the equality will hold if, and only if, f{x) = kg{x), with k real. 
Holder's inequality 

Let f{x) and g{x) be any two functions for which |/(a;)|^ and are integrable 

on [a, b] with p > 1 and + = 1; then, 

''dxf{x)g{x) < ( f\x\f{x)A +( [\x\g{x)A ' . (D.2) 



The equality holds if, and only if, a = /3 , where a and /5 are positive 

constants. 



^The list follows the section on integral inequalities in ref. [Q. The standard reference for 
inequalities from a mathematical viewpoint is B^. 
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Minkowski's inequality 



Let f{x) and g{x) be any two functions for which \ f{x)f and are integrable 

on [a, b] with p > 0; then, 

(^y dx|/(x)+^(a:)|^J <(^y dx\fix)fj + i^j dx\gix)fj . (D.3) 
The equality holds if, and only if, f{x) = kg{x) for some non- negative k. 



Jensen's inequality 

Let f{x) and g{x) be defined for a < x < b such that a < f{x) < (3 and g{x) > 0, 
with g{x) 7^ 0. If (f){u) is a convex function on the interval a < u < (3, i.e. if for 
any two points mi,M2 in 

Ui+U2\ ^ + 0(m2) 



then the following relation between /(x), g{x), and holds: 

XN^^/O^iMfO ^ ^ Iadx^{f{x))g{x) 
Jl^dxg{x) J ~ t^dxg{x) 



Due to the importance of convex functions in statistical mechanics, Jensen's in- 
equality is one of the more widely used integral inequalities. For the use of Jensen's 
inequality in statistical mechanics see the paper by Luttinger, ref. . 
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